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SECTION I 


INTRODUCTION 


The FORTRAN Math Library consists of an extensive assortment of subroutines to aid 
the programmer in performing mathematical and trigonometric operations and functions, 
conversions between data types, bit string operations, logical relations, and other functions. 
The math routines included are for single-(real) and double -precision, complex, integer, 
and logical calculations. 

This library may be loaded in either normal or extended mode and will run in the same 

mode. 

SUBROUTINE DESCRIPTIONS 

The descriptions, in Sections IV and V, of the FORTRAN external and intrinsic functions 
and the compiler support subroutines give the name of the subroutine, its purpose, the 
DAP-16 Mod 2 calling sequence, the FORTRAN calling sequence (where appropriate), the 
method used to compute the result, the data types of the arguments and the result (where 
applicable), error messages generated by the subroutine, if any, and other routines used by 
the subroutines, if any. 

APPENDICES 

There are five appendices to this manual. Appendix A lists the contents of the library 
tapes. There are four magnetic tapes and eight paper tapes available. The first two mag- 
netic tapes are source tapes and contain the sources for the Statistical Library, FORTRAN 
Library, and Fixed Point Math Library. The third magnetic tape contains the objects for 
the software version of the three libraries; the fourth tape contains the objects for the 
hardware version of the three libraries. 


* The list of "Other Routines Used" is given in the order in which they are called. If a 
routine is called more than once, it is listed only once, the first time it is called. 
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The eight object paper tapes contain the FORTRAN Library and are labeled: 


LTCF1 

Tape 1 of b 



LTCF2 

Tape 2 of 6 



LTCF3S 

Tape 3 of 6 

Software Version 

LTCF3H 

Tape 3 of 6 

Hardware 

Version 

LTCF4 

Tape 4 of 6 



LTCF5S 

Tape 5 of 6 

Software Version 

LTCF5H 

Tape 5 of 6 

Hardware 

Version 

LTCF6 

T ape 6 of 6 




The tapes are order dependent, as many of the subroutines call other subroutines which 
appear later in the library. The digits 1 through 6 on the label indicate the loading order. 

Appendix B lists the math routines by argument type. 

Appendix C lists the library subroutines by function. 

Appendix D is an alphabetical list of all the subroutines with their entry points, 
approximate storage required, subroutines referenced and the number of times referenced, 
the library tape on which they are located, and the page in this manual on which they are 
described. 

Appendix E lists the error messages produced by the subroutines and the interpretation 
of these messages. 

SYMBOLS 

The following symbols and letters are used in many of the subroutine descriptions: 

* multiplication 

/ division 

**n raised to the exponential power of n 

C complex 

D double -precision 

I integer 

L logical 

R real 

NAMING CONVENTIONS 

The intrinsic and external functions are named according to the American National 
Standards Institute (ANSI) or the Instrument Society of America (ISA) naming rules. 

The compiler support subroutines are named, for the most part, according to the 
following naming convention: The first letter of the name denotes the operation to be per- 
formed (see the list below). It is followed by a dollar sign having no significance and then 
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by two numbers. The first number (see the list below) represents the operand initially in the 
accumulator (except in load operations) and the second number represents the second operand 
or the type of result. If there is a High-Speed Arithmetic Option version of these subroutines, 
an X is appended to the name. 


Operation 

A - Add 
C — Convert 
D — Divide 
E — Exponential 
H — Hold (store) 
L — Load 
M — Multiply 
N — Negate 
S — Subtract 
Z — Zero (clear) 


Argument Type 

1 — Integer 

2 - Real 

3 — Logical 

5 — Complex 

6 — Double-precision 

8 — Double -precision exponent 


Example s 

A$22 — add two real numbers 

D$52 — divide a complex number by a real number 

E$61 calculate the value of a double -precision number to an 

integer power 

M$22X - multiply two real numbers, using the High-Speed 

Arithmetic Option 


LOADING INFORMATION 

There are two sets of library subroutines, one for installation with the High-Speed 
Arithmetic Option and one for those systems without this hardware option. Each set is 
contained on six rolls of paper tape. Customers who purchase the library in source form 
(on magnetic tape) receive both sets of library subroutines. 

The organization of the library is modular, thus making it possible to load only those 
routines which will be used. This concept of modularity extends to the paper tape. If 
complex or double -precision variables are not used, the first two paper tapes are not required. 

Each paper tape has been assembled via the DAP- 16 Mod 2 assembly language and 
should be loaded by the Series 16 Loader, LDR-APM. Refer to the Series 16 Equipment 
O perators' Manual, Order Number BX48, for information concerning loading object paper 
tapes. 
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SECTION II 

USE OF FORTRAN MATH LIBRARY 


DATA TYPES AND REPRESENTATIONS 

The representation of a negative number in any of the following formats (excluding logical) 
is the TWOs complement of the equivalent positive number. The complement is taken for the 
entire representation, including all subfields. The TWOs complement is taken by reversing all 
bits m the representation (ONEs complement) and adding one to the low-order position, propa- 
gating carries as required. 

Integer 

This is a 16 -bit (single -precision only) word with an implied decimal point after bit 16; 
bit 1 is a sign bit (see Figure 2-1). An integer value may range from -32,768 to +32,767. 

Example: +5 = 0 000 000 000 000 101 = '000005 1 
-5=1 111 111 111 in oil = '177773 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 

Li — i i i i i\i i i » i i i t 1 

'-SIGN BIT 


Figure 2-1. Format of Integer 


Real 

This is a 32 -bit word in the format shown in Figure 2. 2. Bit 1 is the sign bit (0 for 
positive, 1 for negative). Bits 2-9 contain a binary number (N) with a maximum decimal value 
of 255 (377 octal) representing the 8-bit characteristic. This number is "biased" by 128 (200 
octal). The remaining 23 bits represent a binary fraction (F) with a value less than 1. The 
value represented is F*2**(N-128). A number is considered "normalized" when the fraction 
F is at least 1/2 (i. e. , the leading bit is set for a positive number). Within this representation 
the largest representable number in normalized form is just under 2**127, or approximately 
10**(38.5). The smallest number is 2**(-129), or approximately 10**(-38. 5). The 23 magni- 
tude bits give a precision of one part in 2**23, or approximately 6. 9 digits of accuracy. Zero 

The apostrophe before a number indicates octal code. 
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is shown by all zeros in these 23 bits. (Throughout this manual the word "real" is used to 
reference real singie-precision numbers.) 

Example: +5. = 0 100 000 111 010 000, 0 000 000 000 000 000 = '040720, 0 
-5. = 1 011 111 000 110 000, 0 000 000 000 000 000 = '137060, 0 

Double- Precision 

This three-word format is identical to the real number format with the exception of an 
additional 16 magnitude bits (see Figure 2-2). The 39 magnitude bits give a precision of one 
part in 2**39, or approximately 11.7 digits of accuracy. This data type should not be confused 
with hardware double-precision. 

Complex 

This is represented by two real number pairs, each having the format of a real number 
(see Figure 2-2). A real number takes two words of storage; the complex format requires 
four words. The first two words represent the real portion of the complex number, and the 
last two words represent the imaginary portion. 

Lo gical 

A logical value is shown as a word of all zeros for false and a value of one for true. In 
logical operations, any nonzero value is interpreted as true. The complement of a logical 
value changes it from 0 to 1 or 1 to 0. 


FIRST WORD: 
(REAL AND 
DOUBLE- 
PRECISION) 


l 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 



1 1 1 1 1 1 1 

J 

1 1 1 L — 1 

i 

v : 'i 

w 


CHARACTERISTIC 
(EXCESS - 128 NOTATION) 
-SIGN BIT 


MOST SIGNIFICANT 
7 BITS OF THE FRACTION 

L POSITION OF THE BINARY 
POINT 


SECOND WORD: 
(REAL AND 
DOUBLE - 
PRECISION) 


12345 678 9 10 11 12 13 14 15 16 


1 . 

l 1 1 1 1 1 1 

1 i i—L 

_i i 

v 



y 


v 

NEXT MOST SIGNIFICANT 16 BITS OF THE FRACTION 


THIRD WORD: 
(USED FOR 
DOUBLE- 
PRECISION ON LY 


2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 


J L 


l l I I 1 I I I 1 1 L 




LEAST SIGNIFICANT 16 BITS OF THE FRACTION 


Figure 2-2. Format of Real and Double -Precis ion Numbers 
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NORMALIZATION 


A real, double-precision, or complex number is defined as normalized when the fractional 
part has a value between 1/2 and 1. For instance, 3/8 x 2 3 and 3/4 x 2 2 both have the same 
value, but the latter is the normalized form. 

REGISTER USE 

All registers are presumed to be available to the subroutine library, and the user is 
cautioned not to expect any of them to be preserved, whether or not the arguments or results 
are stored m them. That is, any registers not specifically described as containing a particular 
result upon exit from the subroutine must be considered as having become undefined by the 
execution of the subroutine. 

ACCUMULATORS 
Integer Accumulator 

The A-register is used in all integer operations. 

Real Accumulator 

The A- and B-registers are used in all real operations. 

Complex (pseudo) Accumulator 

This four-word area in memory (AC1-AC4) is provided by the library to be used in all 
complex operations. The real portion of the complex number is stored in locations AC1 and 
AC2; the imaginary portion is stored in locations AC3 and AC4. 

Double-Precision (pseudo) Accumulator 

This three-word area in memory (AC1-AC3) is provided by the library to be used in all 
double-precision operations. 

RESULTS 

Results are stored according to their data types. Complex numbers are found in the 
complex accumulator upon exit from any of the compiler support subroutines; double -precision 
numbers are found in the double -precision accumulator; real numbers are found in the A- and 
B-registers; and integer and logical values are found in the A-register. 
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SECTION III 

DAP- 16 MOD 2 PROGRAMMING INFORMATION 


LIBRARY CALLS 

The DAP- 16 Mod 2 calling sequences for entry into the subroutines are shown in the 
descriptions in Sections IV, V and VI. When the FORTRAN compiler encounters either a 
function reference or a call to a subroutine, the following steps are performed: 

1. A call to the function or subroutine is generated. 

2. The address of each argument is determined and saved, in the order in 
which it is retrieved. In the case of expressions, this address is the 
location containing the current value of the expression. 

3. If there are two or more arguments, the final address is followed bv a 
word of zeros to serve as an argument list terminator. 


The code generated by a subprogram definition written in FORTRAN includes a call to 
the special subroutine F$AT (Argument Transfer; refer to Section VI). This call immediately 
follows the entry point and in turn is followed by a word containing a count of the number of 
arguments as defined in the definition statement, followed by that number of words. The 
F$AT subroutine fills in those words with the argument addresses (from the call to the 
subprogram) and sets the return to the word following the argument terminator word (zeros). 

All levels of indirect addressing are removed in passing these addresses. In the case of a 
single argument, the terminator word is eliminated, the argument to F$AT shows a single 
argument, and the search for the terminator is not performed. 

Null arguments may be included in a calling sequence by use of DAC *0 as the address in 
the call. Subroutines serviced by F$AT find the address DAC*0 placed in the list of addresses 
and therefore know that the parameter was null. It is equally effective to use a DAC*PTR, 
where PTR is a DAC*0. This permits a dummy argument to be null, i.e. , an argument 
passed through an intermediate subroutine call. 

The DAP- 16 Mod 2 programmer can generate his own code, performing the same functions 
as the F$AT subroutine. 

Some of the FORTRAN Math Library subroutines have additional arguments in the A- and 
B-registers, or the C-register, or the pseudo -accumulators AC1-AC4. When this is the case, 
the description references an "implicit" argument, i.e., one whose address is not explicitly 
part of the calling sequence. 
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Compiler support subroutines are those which are not normally explicitly called by the 
FORTPAN programmer. For example, the statement 
Z = X + Y 

produces the following DAP- 16 Mod 2 code: 


CALL 

L$22 \ 

loads the value of X in the A- and 

B-registers 

DAC 

X J 


CALL 

A$22 \ 

adds the value of Y to the A- and 

B-registers 

DAC 

Y J 



CALL 

H$22 \ 

stores the result in the A- and B- 

•registers in 

DAC 

z 5 

lo cation Z 



Subroutines L$22, A$22, and H$22 are compiler support subroutines. They may be called 

explicitly by the FORTRAN programmer, if desired, as follows: 

CALL L$22(X) 

CALL A$22(Y) 

CALL H$22(Z) 

To perform the same function as the statement Z = X + Y and to generate the same code. 

Any of the compiler support subroutines may be called by the FORTRAN programmer m 

the following manner: 

CALL ROUTINE NAME (ARG1) 

CALL ROUTINE NAME (ARG1 , ARG2) 

CALL ROUTINE NAME (ARG 1 , ARG2, . . . , ARGn) 

EXAMPLES OF CALLS TO LIBRARY 

CALL M$55 

DAC ARG 1 

Return 

This call enters the complex multiplication subroutine, multiplying the contents of the 

complex pseudo -accumulator by the complex value in locations ARG 1 -ARG 1+ 3 m the standard 

format for complex numbers. The result is stored in the complex accumulator (AC1-AC4), 

and any of the other registers should be presumed to have become undefined. 

CALL AMINO 

DAC I 

DAC J 

DAC K 

OCT 0 

Return 

This subroutine compares the three integer arguments I, J, and K (no implicit arguments) 
and returns with the value of the smallest of these, converted to data type real, in the A- and 
B-registers. Other registers are now nresumed to be undefined. 
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SECTION IV 

INTRINSIC AND EXTERNAL FUNCTIONS AND SUBROUTINES 


This section describes the mathematical and trigonometric functions and special 
FORTRAN subroutines, arranged in alphabetical order by subroutine name. 


4-1 


AM74 



Purpose 


To generate the absolute value of a real number. 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 

RffnlR 


Other Routines 
Used 


CALL ABS 

DAC ARG1 (a real number) 
(Return) 


ABS(R) 


This subroutine checks the real agrument, ARG1, for its algebraic 
sign. If the sign is negative, the TWOs complement of ARG1 is 
calculated. If the sign is positive, the number remains unchanged. 


This absolute value function of a real number gives a real 
result. 


L$22, N$22 
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AIMAG 


Purpose 


DAP Calling 
Sequence 


FORTRAN Refer 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To obtain the imaginary part of a complex argument and convert it 
to real format. 


CALL AIMAG 

DAC ARG1 (a complex number) 
(Return) 


AIMAG (C) 


The complex argument, ARG1, is placed in the complex accumu- 
lator. The imaginary part of the complex number (AC3 and AC4) 
is then loaded into the A- and B-registers. 


The imaginary part of the complex argument, ARG1 , is converted 
to a real number and placed in the A- and B-registers. 


L$55, L$22, AC 3 
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AINT 


Purpose 


To truncate the fractional bits of a real number. 


DAP Calling 
Sequence 


CALL AINT 

DAC ARG 1 (a real number) 
(Return) 


FORTRAN Reference AINT(R) 

Method A constant (2**22) is successively added and subtracted from ARG 1. 

... The ava iiable precision of real numbers is such that the fractional 

part of this result is lost. If ARG1 is negative, its TWOs comple- 
ment is taken before the addition and subtraction take place and it 
is recomplemented before the subroutine exits. The resultant value 
is effectively the largest integer < | ARGl| with the sign of ARG1. 


Data Type of The real argument remains a real number. 

Arguments and 

Results 

Other Routines L$22, N$22, A$22, S$22 

Used 
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ALOG 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To calculate the natural (base e) or common (base 10) logarithm 
of a real number. 


CALL ALOG (or ALOGIO) 
DAC ARG1 (a real number) 
(Return) 


ALOG(R) or ALOGIO(R) 


The log 2 of the argument, ARG1, is computed. This value is 
then converted to the desired base by multiplication by an appro- 
priate constant. 


log ARG 1 = F 1 ;: 

1 

where T = F *F 


*(Cl+T(C3+T(C5+T(C7+T(C9)))))+B-.5 


Cl = 

. 28853901E1 

C3 = 

. 96179665E0 

C5 = 

. 57708664E0 

C 7 = 

. 41 153510EO 

C9 = 

. 34280712E0 



F +v^L 


2 

F is the fractional part of the normalized argument and B is the 
binary exponent of the original argument which has been converted 
to a real number. 


The argument and the results are both real numbers. 


The message "LG" is reported if a negative or zero-valued argu- 
ment is used, and the result is undefined. 


ARG$, C$12, A$22, M$22, S$22, F$ER, H$22, L$22, D$22 
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ALOGX 


Purpose 


To calculate the natural (base e) or common (base 10) logarithm 
of a real number. 


DAP Calling 
Sequence 


CALL ALOGX (or ALOG or ALOG1 0) 
DAC ARG1 (a real number) 
(Return) 


FORTRAN Reference ALOG(R) or ALOGIO(R) 


Method 


log Z = (log Z)*(log 2), where Z=ARG1. Thus for the natural 
A 2 A. 

logarithm, 

In Z = (log Z)*(log 2); for the common logarithm, 

2 2 

log Z = (log., Z)*(logj 0 2). The calculation simplifies in both 
cases to a computation of log 2 Z. Remembering that the floating- 
point number Z can be expressed as Z = F*2**B, where F is the 
fractional part and B the binary exponent of the normalized argu- 
ment Z , 

log 2 Z r (ln(F)/ln (2))+B. 

Now, let F = F*K/K, where K may be the product 
t 

TT K. 

1 

i=l 


such that F*K = 1+G,- where G is positive 
log 2 X = ln(F*K/K) 


+ B 


In (2) 

ln(F*K) - ln(K) 
In (2) 

= ln(F*K) - ln(K) 
ln(2) ln(2) 

= ln(l+G) - ln(K) 


ln(2) 


ln(2) 


+ B then ln(K) is 

• + B 

+ B 


l ln(K ] 


i = l 


G - 1/2G 2 + 1/3G 3 - 


ln(K) 

ln(2) 


+ B 


ln(2) 

Since ln(2) = . 69314718, 

log 2 X = 1.442695141 G - . 7213475704 G 2 + .4808984995 G~ 


ln(K) 

ln(2) 


+ B 
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ALOGX conf 


Data Type of This function with a real argument results in a real number. 

Arguments and 

Results 


Error Messages 


The message "LG" is reported if a negative or zero-valued argu- 
ment is used, and an undefined result is returned. 


Other Routines ARG$, C$12, A$22X, M$22X, S$22X, F$ER 

U sed 
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ALOGIO 


Purpose 


To calculate the common (base 10) logarithm. 


See ALOG or ALOGX . 
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AMAXO 


Purpose 


To find the maximum real value in a list of integers. 


See MAXO. 
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AMAX1 


Purpose 


To find the maximum real value in a list of real arguments. 


See MAXI. 
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AMINO 


Purpose 


To find the minimum real value in a list of integers. 


See MINO. 
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AMIN1 


Purpose 


To find the minimum real value in a list of real arguments. 


See MIN 1 . 
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AMOD 


Purpose 


To compute the remainder resulting from the division of two real 
numbers. 


DAP Calling 
Sequence 


CALL AMOD 

DAC ARG1 (real dividend) 

DAC ARG2 (real divisor) 

OCT 0 (end of arguments flag) 

(Return) 


FORTRAN Reference AMOD(R, R) 


Method This subroutine divides ARG1 by ARG2 by calling D$22. The functior. 

AMOD (ARG 1 , ARG2) is defined as: 

A1 - (A1/A2) * A2, where A 1 = ARG 1 and A2=ARG2 

(A1/A2) is the integer whose magnitude does not exceed the magni- 
tude of A1/A2 and whose sign is the same as that of A1/A2. 


2 a l a Type This function with two real arguments results in a real number for 

Arguments and a remainder. 

Results 


Other Routines L$22, D$22, AINT, M$22, N$22, A$22 

Used 
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ATAN 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Results 


To calculate the principal value of the arctangent (i. e. , 1st or 4th 
quadrant angle) of a real number or to compute and adjust for 
quadrant the arctangent of a real number expressed as a ratio 


(X/Y). 


or CALL AT AN 2 

DAC ARG1 (both arguments 

DAC ARG2 are real numbers 

OCT 0 (end of arguments flag) 
( Return) 


ATAN(R) or ATAN2(R, R) 


CALL ATAN 

DAC ARG1 (a real number) 
(Return) 


For ATAN, let N = ABS(ARGl).. The arctangent of N is evaluated by 
dividing the total range 0<N<10**75 into three intervals: 

If N <10**(-B), ATAN(N) =N 

If N > 1 0**10, ATAN(N) = pi/2 

If 10** (-8) < N < 10**10, 

ATAN(N) = base angle + P(Z) 

= base angle + C1*Z+C3*Z**3+C5*Z**5. . . 

If N <1/2, Z - N and base angle *= 0 

If N<2, Z = (N-1)/(N+1) and base angle = pi/4 

If N<2, Z = (-1/N) and base angle = pi/2 

For ATAN2, the arctangent of the quotient of ARG1/ARG2 (ARG1 = 
side opposite, ARG2 = side adjacent, or sin/cos) is computed as 
in ATAN and adjusted for quadrant by examination of the signs of 
the numerator and denominator. 


Quadrant 

ARG1 

ARG2 

Quotient 

Results (radians 

1 

+ 

+ 

0 to 00 

0 to pi/2 

2 

+ 

- 

- to 0 

pi/2 to pi 

3 

- 

- 

0 to 00 

-pi to -pi/2 

4 

_ 

+ 

- to 0 

-pi/2 to 0 
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ATAN cont 


D M a T YP e ° f This arctangent function of a real number results in a real number 

Arguments and 

Results 


Other Routines ARG$, D$22, N$22, M$22, A$22, S$22 

Used 
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ATAN2 


Purpose 


To calculate the arctangent as the quotient of two real numbers. 


See ATAN. 
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Purpose 


To generate the absolute value of a complex number. 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL CABS 

DAC ARG1 (a complex number) 
(Return) 


CABS(C) 


The argument is squared and its square root is taken to arrive at 
its absolute value; e. g. , if ARG1 = X+IY, 

CABS(ARGl) = SORT (X**2+Y**2). 


This absolute value function of a complex number gives a real 
result. 


F$AT, SUB$, L$22 , M$22, H$22, A$22, SORT 
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Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To calculate the cosine of a complex number with the real part 
in radian measure. 


CALL CCOS 

DAC ARG 1 (a complex number) 
(Return) 


CCOS(C) 


The cosine function is transiormed into the sine function by use of 
the trigonometric identity COS (Z) - SIN (Z + pi/2), where Z=Y+IY. 
SIN (Z+pi/2) is then evaluated. 


This cosine function of a complex number results in a complex numbe 


F$AT, L$55, A$55, H$5B, CSIN 
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Purpose 


To calculate the exponential of a complex number with the imaginary 
part in radian measure. 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL CEXP 

DAC ARG1 (a complex number) 
(Return) 


CEXP(C) 


The following algorithm is used to calculate the value of e**ARGl, 
where ARG1 is a complex number: 

If ARG1 = X+IY, 

e**(X+IY) = (e**X) * (e**IY) = (e**X) *COS(Y) + I % (e**X) *SIN(Y) 


This function raises e to a complex power and gives a complex result. 


F$AT, SUB$, EXP, H$22, COS, M$22, SIN, L$55 
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CLOG 


Purpose 


To calculate a particular value of the natural logarithm (base e) 
of a complex number. 


DAP Calling 
Sequence 


CALL CLOG 

DAC ARG1 (a complex number) 
(Return) 


FORTRAN Reference CLOG(C) 

Method The following algorithm is used to calculate In(ARGl), where 

ARG 1 = X + IY: 

In ^X 4 - IY) - P + T ( A ) 

where R = In (X**2+ Y **2) **. 5 - 1/2 In (X : -'- ;: 2 + Y **2) 

4) - (TAN**- 1)(Y /X) = 4> + 2Kpi 
where K = 0, ± 1, ± 2, ... 

A particular value for <p is chosen such that -pi<(f>< pi by enter- 
ing the arctangent routine ATAN2. 


Data Type of This logarithm function of a complex number gives a complex result. 

Arguments and 

Results 


Other Routines F$AT, L$22, M$22, H$22, A$22, ALOG, ATAN2, L$55, SUB$ 

Used 
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CMPLX 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To combine two real numbers into one complex quantity. 


CALL CMPLX 

DAC ARG1 (a real number) 

DAC ARG2 (a real number) 

OCT 0 (end of arguments flag) 

(Return) 


CMPLX (R, R) 


The first real argument (ARG1) is stored in the real portion of the 
complex accumulator (AC 1 and AC2). The second real argument 
(ARG2) is stored in the complex portion of the complex accumulator 
(AC 3 and AC4). 


The two real arguments are combined into one complex number and 
stored in the complex accumulator. 


F$AT, SUB$, L$22, H$22, L$55 
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CONJG 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To obtain the conjugate of a complex number. 


CALL CONJG 

DAC ARG1 (a complex number) 
(Return) 


CONJG(C) 


This subroutine reverses the sign of the imaginary part of the 
complex argument (ARG1). 


The complex argument in this function remains a complex number. 


F$AT, SUB$, L$22, H$22, N$22, L$55 
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Purpose 


To calculate the cosine of a real number expressed in radians. 


See SIN. 
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CSIN 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To calculate the sine of a complex number with the real part in 
radian measure. 


CALL CSIN 

DAC ARG1 (a complex number) 
(Return) 


CSIN(C) 


The sine function of the complex number ARG1 (X+IY) is computed 
as follows: 

SIN (X+IY) = SIN (X ) * COSH (Y) + I * (COS(X) * SINH(Y) ) 

where SINH(Y) = 1/2 * (E**Y-E**-Y) 

COSH(Y) = 1/2 * (E**Y+E**-Y) 


The argument and the result of this function are complex numbers. 


F$AT, SUB$, EXP, H$22, L$22, D$22, A$22, SIN, M$22, S$22, COS, 
L$55 
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Purpose 


To calculate the square root of a complex number. 


DAP Calling 
Sequence 


FORTRAN Referenc e 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL CSQRT 

DAC ARG1 (a complex number) 
(Return) 


CSQRT(C) 


If the complex argument is positive, (A+B)**. 5 - C+DI is determined 
as follows: 

C = ( ( (A**2+B**2)*®-. 5+A)/2)**. 5 
D = B/(2*C) 

If the argument is negative, ABS(D) = ( ( (A**2+B**2) -A)/2 )**. 5. 

The sign of the real part of the result will be positive and the sign 
of the imaginary part of the result will be the same as the sign of 
the imaginary part of the argument. That is, the results will lie 
in quadrants I or IV of the complex plane. 


This square root function of a complex number results in a complex 
number. 


F$AT, SUB$, CABS, H$22, ABS, A$22, M$22, SORT, L$22, D$22, 
L$55 
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Purpose 


To generate the absolute value of a double -precision number. 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL DABS 

DAC ARG1 (a double -precision number) 
(Return) 


DABS(D) 


This subroutine checks the double -precision argument, ARG1, for 
its algebraic sign. If the sign is negative, the TWOs complement o 
ARG1 is calculated. If the sign is positive, the number remains 
unchanged. 


This function with a double -precision argument results in a 
double -precision number. 


F$AT, L$66, N$66 
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DATAN 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To calculate the arctangent of a double -precision number. 


CALL DATAN 

DAC ARG1 (a double -precision number) 
(Return) 


DATAN (D) 


The principal value is computed. See "Method" for ATAN. 


This function with a double -precision argument results in a double - 
precision number. 


F$AT, DABS, H$66, C$81, L$66, A$66, D$66, M$66, N$66 
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DATAN2 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Re suits 


Error Messages 


Other Routines 
Used 


To calculate the arctangent of the quotient of two double -precision 
numbers. 


CALL DAT AN 2 

DAC ARG1 (a double -precision number (X) ) 

DAC ARG2 (a double -precision number (Y) ) 

OCT 0 (end of arguments flag) 

(Return) 


DATAN2(D,D) 


The arctangent of the quotient (X/Y) is adjusted for the quadrant 
by examining the signs of the numerator and denominator. See 
"Method" for AT AN. 


This arctangent function of a double -precision quantity gives a 
double -precision result. 


The error message "DT" is reported if the second argument is 
zero. The result in the double -precision accumulator is undefined. 


F$AT, L$66, H$66, F$ER, D$66, DATAN, S$66, a$66 
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Purpose 


To convert a real number to double -precision format. 


DAP Calling 
Sequence 


FORTRAN Refe rence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL DBLE 

DAC ARG1 (a real number) 
(Return) 


DBLE(R) 


This subroutine stores the real argument, ARG1, in AC1 and AC2. 
A word of zeros is appended to the real number as the least signifi- 
cant word of the double -precision fraction and stored in AC3. 


The real argument is converted to a double -precision number. 


F$AT, L$22, C$26 
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Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
U sed 


To calculate the cosine of a double -precision number expressed in 
radians. 


CALL DCOS 

DAC ARG1 (a double -precision number) 
(Return) 


DCOS(D) 


The cosine function is transformed into the sine function using the 
trigonometric identity COS (X) = SIN (pi/2+X). SIN (pi/2+X) is then 
evaluated, with X = ARG1. 


This function with a double -precision argument gives a double - 
precision result. 


F$AT, L$66, A$66, H$66, DSIN 
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J urpo£ 


To calculate e**x, where x is a double-precision number. 


DAP Calling 
Sequence 


CALL DEXP 

DAC ARG1 (a double -precision number) 
(Return) 


FORTRAN Reference DEXP(D) 


Method 


Data Type of 
Arguments and 
Results 


In calculating e**ARGl, the following method is used: e**ARGl = 
2**(ARGl*log 2 (e)) = 2**(I+F), where land F are the integer and 
fractional portions, respectively, of the product ARGl*log 2 (e). 


This function raises e to the power of a double -precision argument 
and gives a double-precision result. 


Other Routines 
Used 


F$AT, L$66, M$66, H$66, C$61, 
A$81 


C$16, N$66, A$66, S$66, D$66, 
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DIM 


Purpose 

To compute the 

DAP Calling 

CALL 

DIM 

Sequence 

DAC 

ARG1 

DAC 

ARG2 


OCT 

0 


To compute the positive difference between two real arguments. 


(a real number) 

(a real number) 

(end of arguments flag) 


(Return) 

FORTRAN Reference DIM(R, R) 


Method 


ARG1 - ARG2 is computed. If the result is positive, this value is 
the result given. If ARG1 - ARG2 is a negative quantity, the 
result of this function is zero. 


Data Type of 
Arguments and 
Results 


This routine to calculate the difference between two real numbers 
results in a real number. 


Other Routines L$22, S$22 

Used 
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DINT 


Purpose 

DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To truncate the fractional bits of a double-precision number. 


CALL DINT 

DAC ARG1 (a double -precision number) 
(Return) 


DINT(D) 


A constant (2**38) is successively added and subtracted from the 
argument, ARG1. The available precision of double -precision 
numbers (39 bits) is such that the fractional part of this result is 
lost. If ARG1 is negative, its TWOs complement is taken before 
the addition and subtraction take place and it is recomplemented 
before the subroutine exits. The resultant value is effectively 
the largest integer <|ARG1 | with the sign of ARG1 . 


The double -precision argument after truncation remains a double- 
precision number. 


L$66, N$66, A$66, S$66, AC1 
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DLOG 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To calculate the natural (base e) logarithm of a double - precision 
number. 


CALL DLOG 

DAC ARG1 (a double -precision number) 

(Return) 


DLOG(D) 


This routine is also used by DLOG2 and DLOGIO. Log A (X), where 
X = ARG1, is calculated as log 2 (X)/log 2 (A). To calculate log 2 (X), 
X is considered as the number F *(2**B)^ where 1/2 <F^< 1. Log 2 
(X) = log 2 (F 1 ) + the binary exponent of F 1 , and log 2 (F ) = 1/2 + 
C1*Z + C3(Z**3) + ... where 

Z = (F * ' Cl = 2. 885390081845024D0 

(Fl + V2) C3 = . 96 17966484737 566 DO 
C5 = . 57 7086 6246 3 95 3 5D0 
Cl = .411 53 50984 5700 17 DO 
C9 = . 3428071228932386D0 

This natural logarithm function of a double-precision argument 
results in a double-precision number. 


The message "DL" is reported if a negative or zero -valued argu- 
ment is found. The result in the double -precision accumulator is 
undefined. 


F$AT, DLOG2, M$66 
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Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 

Error Messages 


Other Routines 
Used 


To calculate the common (base 2) logarithm of a double -precision 
number. 


CALL DLOG2 

DAC ARG1 (a double -precision number) 
(Return) 

DLOG2(D) 


This routine is used by DLOG and DLOGIO to calculate log ? (X), 
where X is equal to F**(2**B) and 1/2 sF< 1. See "Method" for 
DLOG. 


This common logarithm function with a double -precision argument 
results in a double-precision number. 


The message "DL" is reported if a negative or zero-valued argu- 
ment is found. The result is undefined. 


F$AT, L$66, F$ER, C$81 , C$16, H$66, Z$80, A$66, S$66, D$66, 
M$66 
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DLOGIO 


Pu rpo s e 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 

Results 


Error Messages 


Other Routines 
Used 


To calculate the common (base 10) logarithm of a double-precision 
number. 


CALL DLOGIO 

DAC ARG1 (a double-precision number) 

(Return) 


DLOGIO(D) 


See "Method" for DLOG. 


This logarithm function with a double -precision argument results 
in a double-precision number. 


The message "DL" is reported if a negative or zero-valued argu- 
ment is found. The result is undefined. 


F$AT, DLOG2, M$66 
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DMAX1 


Purpo s 6 To find the largest value in a list of double -precision arguments. 


DAP Calling 

CALL 

DMAX1 


Sequence 

DAC 

ARG1 

(first double -precision argument) 


DAC 

ARG2 

(a double -precision number) 


DAC 

ARGn 

(last double -precision argument) 


OCT 

0 

(end of arguments flag) 


(Return) 

FORTRAN Reference DMAX 1 (D , D , . . . , D ) 


Method 


Compare the arguments and retain the largest value. 


Data Type of 
Arguments and 
Results 


The largest double -precision argument is stored in the double - 
precision accumulator. 


Other Routines L$66, H$66, S$66 

Used 
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DMIN1 


Purpose 

DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To find the smallest value in a list of double -precision arguments. 


CALL 

DMIN1 


DAC 

ARG1 

(a double-precision argument) 

DAC 

ARG2 

(a double -precision argument) 

DAC 

ARGn 

(last double-precision argument) 

OCT 

0 

(end of arguments flag) 

(Return) 


DMIN1 

(D,D, . . 

• . D) 


Compare the arguments and retain the smallest value. 


Both of the arguments are double-precision and the result of this 
function is a double-precision number. 


L$66, H$66, S$66 
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DMOD 


Pu rpo s e 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To compute the remainder resulting from the division of two double - 
precision numbers. 


CALL DMOD 

DAC ARG1 (a double pr ecision number) 

DAC ARG2 (a double -p recision number) 

OCT 0 (end of arguments flag) 

(Return) 


DMOD(D, D) 


This subroutine divides ARG1 by ARG2 by calling D$66. The 
function DMOD (A1.A2) is defined as A1 -(A 1/A2)*A2, where (A1.A2) 
is the integer whose magnitude does not exceed the magnitude of 
A1/A2 and whose sign is the same as that of A1/A2. 


This function with two double-precision arguments results in a 
double-precision number for a remainder. 


F$AT , L$66, D$66, H$66, DINT, M$66, S$66, N$66 
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DSIGN 


Purpose To generate a value consisting of the sign of the second double- 

_ precision argument and the magnitude of the first double-precision 

argument. 

CALL DSIGN 

DAC ARG1 (a double -precision number) 

DAC ARG2 (a double -precision number) 

OCT 0 (end of arguments flag) 

(Return) 


FORTRAN Reference DSIGN(D, D) 

Method ARG2 is tested for its algebraic sign and, depending on the sign of 

_ ARG1, the procedure is as follows: 

ARG 1 ARG2 Result 

+ + j ARG 1 j 

- | ARG1 | 

+ + + | ARG1 | 

+ - - | ARG 1 | 

Data Type of Both arguments for this call are double-precision numbers and the 

Arguments and result is a double -precision number. 

Results 


Other Routines F$AT, L$66, N$66 

Used ” 


DAP Calling 
Sequence 
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DSIN 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To calculate the sine of a double -precision number expressed in 
radians . 


CALL DSIN 

DAC ARG1 (a double -precision number) 
(Return) 


DSIN (D) 


An arbitrary angle X expressed in radian measure can be reduced 
to the range 0 < Y < Jp- through the relation X = Y + N(pi/2). 
Adjustment is made for quadrant before using a modified Taylor's 
expansion. 


This sine function with a double -precision argument results in a 
double -precision number. 


F$AT, DABS, M$66, H$66, C$61, C$16, N$66, A$66, MOD, 
L$66, S$66 
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Purpose 


To calculate the square root of a double -precision number. 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 


CALL DSQRT 

DAC ARG 1 (a double -precision number ) 

(Return) 


DSQRT (D) 


A first approximation to the double -precision square root of the 
double -precision argument is obtained by calling the real square 
root routine (SQRT ). One more Newton-Raphson iteration is the 
made to achieve full double -precision accuracy. 


This square root function of a double -precis ion argument result 
in a double -precision number. 


F$AT , L$66, C$62, H$22, SQRT, C$26, H$66, D$66, A$66, 
A$81 
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Purpose 


To calculate e*-'x, where x is a real number. 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Error Message 


Other Routines 
Used 


CALL EXP 

DAC ARG 1 (a real number) 
(Return) 


EXP(R) 


e**ARGl - 2-*(ARGl * log^e)) = 2 * * (1+ F ) , where I is the integer 
and F is the fractional portion of the product ARG 1 log (e). The 
value of F is used to define the quantities I 1 , F(l), and 2 F(2): 


F 

r 

F(l) 

F (2 ) 

-1 < F < -1/2 

I - I 

1/4 

F + 3/4 

-1/2 < F < 0 

I - I 

3/4 

F + 1 /4 

0 < F < 1/2 

I 

1/4 

F - 1/4 

1/2 <F < 1 

I 

3/4 

F - 3/4 

F rom the above 

table, e* 

-ARG 1=2- 

-(F+F1+F2; 


where 

2*»F2 = e**(F2*ln(2)) = e**F = (A(F) )/(A(F)-B(F) ) 
A(F) = C 1 + (F*F), B(F) = C2*F 


R 

This exponential function with a real argument (e ) results in a real 
number. 


When overflow occurs, the error message "EX" is reported and the 
answer returned is the maximum value possible (1. 7E38). When 
underflow occurs, the value 0 is returned without an error message. 


ARG$, N$22, M$22, S$22, A$22, D$22, F$ER 
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Purpose 


To convert an integer argument to real format. 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL FLOAT 

DAC ARG1 (an integer value) 
(Return) 


FLOAT (I) 


This routine extracts the integer and converts it to real format, 
leaving the result in the A- and B-registers. 


This routine converts an integer argument to a real number. 


C$12 
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IABS 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


To generate the absolute value of an integer. 


CALL IABS 

DAC ARG1 (an integer value) 
(Return) 


IABS (I) 


This subroutine checks the integer argument, ARG1 , for its algebraic 
sign. If the sign is negative, the TWOs complement of ARG1 is cal- 
culated. If the sign is positive, the number remains unchanged. 


This absolute value function with an integer argument results in an 
integer. 
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IDIM 


Purpose 


To compute the positive difference between two integer arguments. 


DAP Calling 
Sequence 


CALL IDIM 

DAC ARG1 (an integer value) 

DAC ARG2 (an integer value) 

OCT 0 (end of arguments flag) 

(Return) 


FORTRAN Reference IDIM(I, I) 


Method 


Compute DIF = ARG1-ARG2. If DIF is positive, the result of this 
function is the value of DIF. If DIF is negative, the result of this 
function is zero. 

DIF = ARG 1 - MIN (ARG 1 , ARG2 ) 


Data Type of The result of this function with two integer arguments is an integer. 

Arguments and 

Results 
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IDINT 


Pprppse To truncate the fractional bits from a double -precision argument, 

thus converting it to integer format. 

See IF IX . 
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IFETCH 


Purpose 


To fetch the contents of the memory location specified by ARG1. 


DAP Calling 
Sequence 


CALL IFETCH 
DAC ARG 1 
(Return) 


FORTRAN Reference IFETCH(ARGl) 


Method The A-register is loaded with the contents of the location specified 

by ARG 1 . 


Other Routines ARG$ 

Used ~~ 
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IFIX 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To truncate the fractional bits from a real or double -precision 
argument, thus converting it to integer format. 


CALL IFIX (or CALL INT ) 

DAC ARG1 (a real number) 

(Return) 


or 

CALL IDINT 

DAC ARG1 (a double -precision number) 
(Return) 


IFIX (R), INT(R), IDINT (D) 


This subroutine truncates the fractional bits of ARG1, shifts it to 
the right until the binary point is at the end of the register, and 
normalizes the result. It then uses the characteristic to scale 
the value to an integer. 


If either IFIX or INT is called, the argument is a real number and 
the result is an integer. If IDINT is called, the argument is a double- 
precision number and the result is an integer. 


L$22 , C$21 
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INT 


Purpose 


To truncate the fractional bits from a real argument, thus convert- 
ing it to integer format. 


See IFIX . 
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ISIGN 


Purpose 


To generate a value consisting of the sign of the second integer 
argument and the magnitude of the first integer argument. 


DAP Calling 
Sequence 


CALL ISIGN 

DAC ARG1 (an integer value) 

DAC ARG2 (an integer value) 

OCT 0 (end of arguments flag) 

(Return) 


FORTRAN Reference ISIGN(I, I) 


Method 


ARG2 is tested for its algebraic sign and, depending on the sign of 
ARG1, the procedure is as follows: 

ARG 1 ARG2 Result 

+ + + IARG1 I 

+ - - |ARG1 | 

+ + IARG1 | 

- | ARG1 | 


Data Type of Both arguments and the result are integers. 

Arguments and 

Results 
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Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Other Routines 
Used 


To store the contents of the second argument in the location speci- 
fied as the first argument. 


CALL ISTORE 

DAC ARG 1 (target word address) 

DAC ARG 2 (word to be stored) 

OCT 0 (end of arguments flag) 

(Return) 


ISTORE (ARG 1 , ARG2) 


Fetch the target word address (ARG1) and save it. 

Fetch the word to be stored (ARG 2) and use it to replace the 
contents of the target location. Effectively, the contents of ARG2 
are stored in location ARG1. 


F$AT 
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LOC 


Purpose 


To determine the address of the argument. 


DAP Calling 
Sequence 


CALL LOC 
DAC ARG1 
(Return) 


FORTRAN Reference LOC (ARG 1 ) 


Method 


Fetch the argument address (direct or indirect) and load it into 
the A-register. 
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MAXO 


Purpose 

To find the largest 


with this 

value or c 

DAP Calling 

CALL 

MAX 0 

Sequence 

DAC 

ARC 1 

DAC 

ARC 2 


DAC 

ARGn 


OCT 

0 


(Return 

) 

FORTRAN Reference 

MAXO (I, 

I, ... ,D 


value in a list of integer arguments and exit 
onve rt it to real format (AMAXO) and exit. 


(or AMAXO) 

(integer value) 

(intege r value ) 

(last integer argument) 
(end of arguments flag) 


or AMAXO (I, I, .... I) 


Method This subroutine compares the arguments and retains the largest 

value. If AMAXO is called, the result is converted to real by- 
calling FLOAT before the subroutine exits. 


Data Type of 
Arguments and 
Re suits 


The arguments are integers in either call (MAXO or AMAXO). The 
result is integer if MAXO is called; the result is a real number if 
AMAXO is called. 


Othe r Routines FLOAT 

Used 


4-54 


AM74 



MAXI 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To find the largest value in a list of real arguments and exit with 
this value or convert it to an integer (MAXI) and exit. 


CALL 

MAXI 

(or AMAX1) 

DAC 

ARG 1 

(a real number) 

DAC 

ARG2 

(a real number) 

DAC 

ARGn 

(last real argument) 

OCT 

0 

(end of arguments flag) 


(Return) 


MAX1(R,R R) or AMAX1(E,R, . . . , R) 


This subroutine compares the arguments and retains the largest 
value. If MAXI is called, the result is converted to integer by 
calling IFIX before the subroutine exits. 


The arguments are real numbers in either call (AMAX1 or MAXI). 
The result is real if AMAX1 is called; the result is an integer if 
MAXI is called. 


L$22, H$22, S$22, IFIX 
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MINO 


Purpose To find the smallest value in a given set of integers and exit with 

this value or convert this value to a real number and exit. 


DAP Calling 

CALL 

MINO 

(or AMINO) 

Sequence 

DAC 

ARC 1 

(an integer value) 

DAC 

ARG2 

(an integer value) 


DAC 

ARGn 

(last integer argument) 


OCT 

0 

(end of arguments flag) 


(Return) 

FORTRAN Reference MINO (I, I, . . . , I) or AMINO (I, I, . . . , I) 


Method This subroutine compares the arguments and retains the smallest 

value. If AMINO is called, the result is converted to a real number 
before the subroutine exits. 


Data Types of 
Arguments and 
Results 


The arguments are integers in either call (MINO or AMINO). The 
result is integer if MINO is called; the result is a real number if 
AMINO is called. 


Other Routines FLOAT 

Used 
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MINI 


P urpose To find the smallest value in a list of real arguments and exit with 

this value (AMIN 1 ) or convert it to an integer (MINI) and exit. 


DAP Calling 

CALL 

MINI 

(or AMIN 1 ) 

Sequence 

DAC 

ARG1 

(a real number) 


DAC 

ARG2 

(a real number) 


DAC 

ARGn 

(last real argument) 


OCT 

0 

(end of arguments flag) 


(Return) 


FORTRAN Reference MINI (R, R, . . . , R) or AMIN1 (R, R, . . . , R) 


Method 


Compare the arguments and retain the smallest value. 


Data Type of The arguments are real numbers for either call (MINI or AMIN 1 ) 

Arguments and The result is real if AMIN1 is called; the result is integer if MINI 

Results is called. 


Other Routines L$22, H$22, S$22, IFIX 

Used 
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MOD 


Purpose 

To compute the remainder resulting from the division of two integers. 

DAP Calling 

CALL MOD 

Sequence 

DAC ARG 1 (an integer value) 

DAC ARG 2 (an integer value) 

OCT 0 (end of arguments flag) 

(Return) 

FORTRAN Reference 

MOD (I, I) 

Method 

This subroutine divides ARG 1 by ARG2 by calling D$ll. The function 
MOD(Al,A2) is defined as A1 - (A1 /A2 )*A2, where (A1/A2) is the 
integer whose magnitude does not exceed the magnitude of A1/A2 
and whose sign is the same as that of A1/A2. 

Data Type of 

This function with two integer arguments results in an integer for a 

Arguments and 

Results 

remainder. 

Other Routines 

Used 

D$1 1 , M$1 1 
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Purpose 


To check for an error condition. 


DAP Calling 
Sequence 


FORT RAN Reference 


Method 


Other Routines 
Used 


CALL OVERFL 

D-AC J (an integer value) 

(Return) 


OVERFL(J) 


This subroutine checks error flag AC 5 for a nonzero value, which 
indicates that an entry to the error subroutine, F$ER, was made 
since the last call to OVERFL. If AC5 is nonzero, the variable J 
is set to 1 and AC5 is cleared. If AC5 is zero, J is set to 2. 


AC 5 
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Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To load the real portion of a complex number into the A- and B- 
register . 

CALL REAL 

DAC ARG 1 (a complex number) 

(Return) 

REAL(C) 


This subroutine calls ARG$ to place the complex argument, ARG 1 , 
into the index register. The real portion, i.e. , the first two 
words, of the complex argument is then loaded into the A- and B- 

registers . 


This function of a complex number results in a real number. 


ARG$ 
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SIGN 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To generate a value consisting of the sign of the second real argu- 
ment and the magnitude of the first real argument. 


CALL SIGN 
DAC ARG1 
DAC ARG2 
OCT 0 
(Return) 


(a real number) 

(a real number) 

(end of arguments flag) 


SIGN (R, R) 


ARG2 is tested for its algebraic sign and, depending on the sign of 
ARG1, the procedure is as follows: 

ARG I ARG2 Result 

+ + +| ARG1 | 

+ - - | ARG1 | 

+ + | ARG1 I 

| ARG 1 | 


Both arguments are real numbers and the result is a real number. 


L$22, N$22 
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SIN 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To calculate the sine or cosine of a real number expressed in radians. 


CALL SIN (or COS) 

DAC ARG1 (a real number) 

(Return) 


SIN(R) or COS (R) 


The angle is reduced to the first quadrant by the use of the relation 
X = Y+N*(pi/2) and the identities SIN(Y) = COS(pi/2-Y) and COS(Y) 
= SIN(pi/2-Y). A modified Taylor's expansion is then used to cal- 
culate the sine of the first quadrant angle. 

The cosine function is transformed into the sine function by the use 
of the Identity COS(X) = SIN(pi/2-X); SIN(pi/2-X) is then evaluated, 
wnere X =ARG1 


This sine function with a rea.' argument results in a real number. 


ARG$, N$22, M$22, S$22, A$22 
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Purpose 


To set or reset the pseudo sense lights and switches. 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Other Routines 
U sed 


CALL SLITE 

DAC ARG 1 

(Return) 


(where ARG1 is the address of the variable con- 
taining the sense light number). 


CALL SLITET 
DAC ARG1 

DAC ARG2 

OCT 0 

(Return) 


(or CALL SSWTCH) 

(where ARG1 is the address of the variable con- 
taining the sense light or switch (SSWTCH) 
number to be interrogated, and ARG 2 is the 
address of the location in which to store the 
"set or reset" indicator; (l=set, 2= reset). 


CALL SLITE (I), CALL SLITET(I, J), CALL SSWTCH(I, J) 


SLITE — The ARG$ routine is used to place the variable address in 
the index register. The argument (I) is tested for zero. If zero, all 
sense light positions are reset; otherwise, the sense light specified 
is shifted to its appropriate position and INCLUSIVELY ORed with 
current settings, leaving them undisturbed. 

gLITET -The ARG$ routine is used to place the sense light number 
in the A-register and the location of the variable in the index register. 
If the sense light number is 0, a 2 is inserted into the variable J, 
signifying a reset condition. Otherwise, the sense light bit is moved 
to its proper position in the A-register. A logical AND is executed 
with the sense light register. If the result of the AND is zero, the 
sense light is reset and a 2 is placed in J. If the result of the AND 
is not zero, an EXCLUSIVE OR is carried out with the sense light 
register, resetting the sense light specified and storing a 1 in J to 
signify that the sense light was set on entry. 

gS WTCH - The ARG$ routine is used to place the sense switch number 
in the A-register and the variable location in the index register. If 
the sense switch number is 0 (no real switch), J is set to 1 . If the 
sense switch number is valid (1 to 4), J is set to 1 if the external 
switch is set and set to 2 if the external switch is not set. 


ARG$, L$33 
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To set or reset the pseudo sense lights and switches. 



See SLITE. 
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Purpose 


DAP Calling 
Sequence 


FORTRAN Refe rence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To calculate the square root of a real number, 
a high-speed version, SQRTX. ) 



(This subroutine has 


CALL SORT 

DAC ARG1 (a real number) 

(Return) 


SQRT(R) 


Given the argument N = F(2**e), the mantissa is adjusted so that 
e is even and 1/4 < e < 1. An initial approximation to the square 
root (Y) is chosen as follows: 

Y = 7/8(F) + 9/32 if e <1/2 

Y = 9/16(F) + 7/16 if e > 1/2 

Two Newton-Raphson iterations are then made to obtain full single- 
precision accuracy. 


This square root function of a real number results in a real number. 


The error message "SQ" is reported if a negative argument is found. 
An undefined result is returned in the A-and B-registers. 


ARG$, DIV$, D$22, A$22, F$ER 
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SQRTX 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To calculate the square root of a real number. (This routine re- 
quires the High-Speed Arithmetic Option. ) 


CALL SQRTX (or SQRT) 

DAC ARG1 (a real number) 

(Return) 


SQRT (R) 


Given the argument N=F*(2~*e), the mantissa is adjusted so that 
e is even and 1/4 < e < 1. An initial approximation to the square 
root of ARG1 is chosen as follows: 

ARG 1 = 7 /8(F) + 9/32 if e <1/2 

ARG1 = 9/ 16(F) + 7/16 if e > 1/2 

Two Newton-Raphson iterations are then made to obtain full single- 
precision accuracy. 


This square root function of a real number results in a real number. 


The error message "SQ" is reported if a negative argument is found. 
An undefined result is returned in the A-and B-registers. 


ARG$ , D$22X, A$22X, F$ER 
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Purpose 


To set or reset the pseudo sense switches. 


See SLITE. 
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TANH 


Purpose 


DAP Calling 
Sequence 


FORTRAN Reference 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
U sed 


To calculate the hyperbolic tangent of a real number. 


CALL TANH 

DAC ARG1 (a real number) 

(Return) 


TANH(R) 


TANH = ( e * * (2 *X ) - 1 ) / ( e * *f2 *X )+ 1 ) , where X = ARG 1 . 


This tangent function with a real argument results in a real number. 


L$22, EXP, A$22, H$22, D$22 
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SECTION V 

COMPILER SUPPORT SUBROUTINES 


This section describes the compiler support subroutines, i. e. , those subroutines which 
are not normally explicitly called by the FORTRAN programmer. These subroutines per- 
form conversions between data types, logical relationals, arithmetic operations, and miscel- 
laneous functions. 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To add or subtract real numbers. (This subroutine has a high- 
speed version, A$22X. ) 

CALL A$22 (or S$22) 

DAC ARG2 (a real number) 

(Return) 

A$22 (Add) - The contents of ARG2 are added to the contents of the 
A- and B- registers after both numbers are unpacked and scaled. 
The result is normalized and the characteristic is adjusted. 

S$22 (Subtract ) - The value contained in ARG2 is negated and the 
add routine, A$22, is entered. 


< implicit real argument> ± < real argument> — *- < real result> 


The error message "SA" is reported if an arithmetic overflow 
occurs, i. e. , the result is > 2**127. An undefined result is 
returned. 


ARG$, N$22, F$ER 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To add or subtract real numbers. (This routine requires the High- 
Speed Arithmetic Option.) 


CALL A$22X (A$22, S$22 or S$22X) 

DAC ARG2 (a real number) 

(Return) 


A$22 (Add ) - The contents of ARG2 are added to the contents of the 
A- and B- registers after both numbers are unpacked and scaled. 
The result is normalized and the characteristic is adjusted. 

S$22 (Subtract ) - The value contained in ARG2 is negated and the 
add routine, A$22, is entered. 


< implicit real argument> ± < real argument> -► < real result> 


The error message "SA" is reported if an arithmetic overflow 
occurs, i. e. , the result is > 2**127. An undefined result is 
returned. 


N$22, F$ER 
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Purpose 


To add a real argument to a complex number. 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL A$52 

DAC ARG2 (a real number) 
(Return) 


The following is the algorithm used to compute the operation of 
adding a real argument (ARG2) to the contents of the complex 
accumulator (Y): 

Y + ARG2 = A + B .* I + ARG2 

= (A + ARG2) + B * I) 
where Y = A + B * I 


< implicit complex argument? + < real argument? — ► < complex 
re suit > 


F$AT, H$55, L$22, A$22, H$22, L$55 
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Purpose 


To add complex numbers. 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL A$55 

DAC ARG2 (a complex number) 
(Return) 


The following is the algorithm used in the addition of two complex 
numbers (the contents of ARG2 and the complex accumulator): 

X+ARG2 = (A+B*I) + (M+N*I) = (A+M) + (B+N) * I 

where X = A+B*I and ARG2 = M+N-I 


< implicit complex argument> + < complex argument> — 

< complex re suit > 


F$AT, H$55, SUB$, L$22, A$22, H$22, L$55 
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Purpose 


To add a real number to a double-precision number. 


DAP Calling 
Sequence 


CALL A$62 

DAC ARG2 (a real number) 
(Return) 


Method This subroutine calls DBLE to convert the real argument to a 

double-precision number and calls A$66 to perform the double- 
precision addition. 


Data Type of < implicit double -precision argument> + < real argument> 

Arguments and < double -precision result> 

Results 


Other Routines F$AT, H$66, DBLE, A$66 

Used 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


To add, subtract, multiply, or divide normalized, double- 
precision numbers. (This subroutine has a high-speed version, 
A$66X. ) 

CALL A$66 (or S$66, M$66, or D$66) 

DAC ARG2 (a double-precision number) 

(Return) 


The contents of ARG2 are added to, subtracted from, multiplied 
by, or divided into the contents of the double -precision accumu- 
lator (X). 


Add (A$66) — The numbers are unpacked and scaled to coincident 
places. The addition process takes place (X+ARG2), and the result 
is normalized. 


Subtract (S$66) - The numbers are unpacked and scaled to coincident 
places. The subtraction process takes place (X-ARG2), and the 
result is normalized. 


Multiply (M$66) - X*ARG2 = (X*2**E1) * (Y*2**E2) 

= X*ARG2*2** (E1+E2) 

Let X = (A+B*2** (-N)) 
and ARG2 = (C+D*2** (-N)) 

X*ARG2 = A*C + ((A*D+B*C) - 2**(-N)) 

The term B*D*2** (-2N) is ignored. 

The least significant bits of the product are: 
L*(A*C)+H*(A*D)+H*(B*C) 


Divide (D$66) — The quotient X/ ARG2 is obtained by the binomial 
expansion of 1/X = X ! ' :> ' : (= 1). The high-order and low-order parts 
(H and L) of the quotient are computed as follows: 

(A+B*2**(-N))/(C+D*2**(-N)) = (A+B-A*D/C)/C 

H = ( A+B- A*D/ C) / C 
L = remainder (H)/C 


< implicit double -precision arguments* 


< double -precision 


argument> — •- < double -precision result> 


The error message "AD" is printed if an addition or 
subtraction over /underflow occurs. 

The error message "PZ" is printed if a division by zero 
is attempted. 

The error message "MD" is printed if a multiplication or 
division over /underflow occurs. 
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A$66 cont 


After an error message is reported, the double-precision accumu- 
lator is loaded with the maximum ((2**128)- 1) or minimum 
(2**(- 128)) value (as determined by the correct sign) before 
returning to the calling program. 


Other Routines N$66, F$ER, H$66, L$66, ARG$, AC1, AC2, AC3 

Used 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 

Other Routines 
U sed 


To add, subtract, multiply, or divide normalized, double -prec ision 
numbers. (This routine requires the High-Speed Arithmetic Option.) 


CALL A$66X (or A$66, S$66, S$66X, M$66, M$66X, D$66, D$66X) 
DAC ARG2 (a double -precision number) 

(R eturn) 


The contents of ARG2 are added to, subtracted from, multiplied by, 
or divided into the contents of the double -prec i sion accumulator. See 
A$66, described on the preceding pages, for a detailed description of 
the methods used. 


implicit double -precision argument > 


... ■< double -precision 


argument? — ► < double -precision result> 


See Error Messages for A$66. 


N$66, F$ER, H$66, L$66, ARG$, AC1, AC2, AC3 
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A$81 


Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Re suits 


Error Messages 


Other Routines 
U sed 


To add an integer value (I) to the characteristic of the variable in the 
double-precision accumulator (effectively, multiplication by 2 1 ). 


CALL A$8 1 

DAC ARG2 (an integer value) 
(Return) 


The characteristic (base 2) of the value in the double-precision accum- 
ulator is increased (or decreased) by an integral value, ARG2. For 
example, if ARG2 = 2 and the value in the double -precision accumu- 
lator is 8.0 (2 3,0 ), the result of this call would be 2 J * or 2 ‘ = 

32.0 (8.0*2 2 ). If the absolute value of the result is less than 2**(-128), 
a value of zero is returned. 


< implicit double -precision argument> * (2** 4 integer argument>) — <► 

< double -precision result> 


If there is exponent overflow, an "EQ" error message is reported and 
external locations AC1 and AC2 are loaded with the maximum value 
po s s ible' ((2**128)- 1 ) with the sign of ARG2. 

N$22, F$ER, AC 1, AC2 
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AC1 


(AC2, AC3, AC4, AC5) 


Purpose To assign locations to be used as a double-precision or complex 

accumulator by the FORTRAN library routines. 

JJ se AC1, AC2, AC3: double-preci sion accumulator. 

AC1, AC2: complex accumulator, real portion. 

AC3, AC4: complex accumulator, imaginary portion. 

AC5: error flag. 
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Purpose 


To convert the indirect address of an argument to its corresponding 
direct address. 


DAP Calling 
Sequence 


CALL ARG$ 

DAC* ARG2 (usually a subroutine entry) 
(Return ) 


Method The address of the argument is returned in the index register. This 

subroutine may be used upon entering a subroutine to set up the return 
addre ss. 
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Purpose 

DAP Calling 
Sequence 

Method 


Data Type of 
Arguments and 
Re suits 


Other Routines 
U sed 


To convert an integer to a real number. 

CALL C$12 
(Retu rn ) 


The integer value in the A-register is placed in the B-register and 
the A-register is set to 045600 (octal), representing a characteristic 
such that the number fits the description given for a real number ex- 
cept that it is not "normalized. " A$22 (with argument = 0 
(040000,000000), also unnormalized) is called to normalize the 
re suit. 


The integer value in the A-register is converted to a real number and 
placed in the A- and B-registers. 


A$22, N$22 
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Purpose 


To convert an integer to a double -pre c ision number. 


DAP Calling 
Sequence 

Method 


Data Type of 
Arguments and 
Re suits 


Other Routines 
U sed 


CALL C$16 
(Return) 


The integer in the A-register is normalized and converted to real by 
calling C$12. This real value is then converted to a double-precision 
number by calling C$26. The result is placed in the double-precision 
accumulator. AC 1 contains the contents of the B-register (the real 
exponent), AC 2 contains the contents of the A-register (the most 
significant word of the fraction), and AC3 contains a word of zeros. 


The integer value in the A-register is converted to a double-precision 
number and placed in the double -precision accumulator. 


C$12, C$26 
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Purpose 


To convert a real number to an integer. 


DAP Calling 
Sequence 

Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


CALL C$21 
(Return) 


This subroutine scales the real number in the A- and B-registers to 
23 bits by adding the octal value 045700 (2**22) to truncate the 
fractional part of the real number. The result is in the A-register. 


The real number in the A- and B-registers is converted to an integer 
and returned in the A-register. 


The message "RI" is reported if the integer (I) is too large when con- 
verted from real to integer. The integer must be in the following 
range: -2^ < I < 2^-1. An undefined result is returned in the A- 
register . 


N$22, A$22, F$ER 
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Pur po se 


To convert a real number to a complex number. 


DAP Calling 
Sequence 


CALL C$25 
(Return! 


Method The A- and B- registers are stored in AC1 and AC2, respectively 

(the real part of the complex number), and AC3 and AC4 (the imagin- 
ary part of the complex number) are set to zeros. 


Data Type of 
Arguments and 
Results 


The real argument in the A- and B- registers is converted to a com- 
plex number and stored in the complex accumulator (AC1, AC2, AC3, 
and AC 4). 


Other Routines H$22, CMPLX 

U sed 
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Purpose 


To convert a real number to a double -precis ion number. 


DAP Calling 
Sequence 

Method 


Data Type of 
Arguments and 
Results 


Other Routines 
U sed 


CALL C$26 
(Return) 


The number in the A- and B-registers is placed in AC1 and AC2. 
AC3 is cleared and the routine exits. 


The real number in the A- and B-registers is converted to double- 
precision and placed in the double-precision accumulator. 


AC1, AC2, AC3 
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Purpose 


To convert a double-precision number to an integer. 


DAP Calling 
Sequence 

Method 


Data Type of 
Arguments and 
Re suits 

Other Routines 
Used 


CALL C$61 
(Return) 


This subroutine calls C$62 to convert the number in the double- 
precision accumulator to real and calls C$21 to convert the real 
number to integer. 


The double-precision value in the double-precision accumulator is 
converted to an integer and placed in the A- register. 


C$62, C$21 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To convert a double-precision number to a real number. 


CALL C$62 or CALL SNGL 

(Return) DAC ARG1 (a double-precision 

(Return) number) 


AC 1 and AC2 (the exponent and the most significant part of the 
fraction of the number in the double-precision accumulator) or the 
first two words of ARG 1 (if SNGL is called) are loaded into the A- 
and B- registers. The least significant part of the fraction (AC3. 
word 3) is not considered in the result. 


The double-precision value in the double-precision accumulator or 
in ARG1 is converted to a real number and placed in the A- and B- 
register s. 


L$22, N$66, N$22, L$66, AC1, AC2 
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Pur po se 


DAP Calling 
Sequence 

Method 


Data Type of 
Arguments and 
Results 

Other Routines 
Used 


To convert the exponent of the value in the double-precision accumu- 
lator to an integer. 


CALL C$81 
(Return) 


Extract the characteristic (base 2) from the value in the double- 
precision accumulator (AC1) and convert it to an integer. 


The characteristic of the double-precision argument is converted to 
an integer. 


AC 1 
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D$11 


Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
U sed 


To divide two integers. (This subroutine has a high-speed version, 
D$1 IX. ) 


CALL D$1 1 

DAC ARG2 (integer divisor) 
(Return) 


The numerator (an integer value) should be in the A-register upon 
entrance to this subroutine. If the denominator, ARG2, is zero, 
an overflow occurs and an error message is reported. If both 
arguments are nonzero, the numerator is positioned in the A- and 
B-registers and the division is performed. The results are ex- 
amined for the special case (-32, 768/ -1) which is treated as an over- 
flow. If the re suits are in the range of - 32, 768 to + 32, 767, D$ll 
returns to the calling program with the quotient in the A-register 
and the remainder in the B-register. The integer answer is in the 
A-register. 


< implicit integer argument > / < integer argument> < integer 
re suit > 


The error message "IZ" is reported if a division by zero is at- , 
tempted. The maximum value is output (-32,768 if negative or 
+ 32, 767 if positive) . A division of -32, 768 by -1 also causes "IZ" 
to be reported; D$ll returns a value of + 32, 767, the maximum 
value possible. 


ARG$, F$ER 
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D$11X 


Purpose 


To divide two integers. (This routine requires the High-Speed 
Arithmetic Option. ) 


DAP Calling 
Sequence 


CALL D$1 1 X (or D$ll) 

DAC ARG2 (integer divisor) 

(Return) 


Method 


See "Method" for D$ll. 


Data Type of < implicit integer argument> / < integer argument> — < integer 

Arguments and re suit > 

Results 


Error Messages 


See "Error Messages" for D$ll. 


Other Routines ARG$, F$ER 

U sed 
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Purpose 


To divide two real numbers. (This subroutine has a high- 
speed version, D$22X. ) 


See M$22 . 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To divide two real numbers. (This subroutine requires the High- 
Speed Arithmetic Option. ) 


CALL D$22X (or D$22) 

DAC ARG2 (the real divisor) 

(Return) 


This subroutine divides the real number in the A- and B-registers 
(X) by the real argument, ARG2 (Y). The division is performed by 
multiplying X by the reciprocal of Y, i.e., X*l/Y. Newton's 
method for 1/Y is: 

R( 1) = R ( 0 ) * (2-R(0)*Y 

where 

R(0) - i/H(Y), H(Y) being the high-order 15 bits of Y 
X * (1/Y) = X * R( 1 ) - X * R(0) * (2 -R(0)*Y) 


< implicit real argument> / < real argument> < real result> 


A !! DZ" error message is typed if division by zero is attempted. 
A value of 0 is returned if the dividend is also 0. The signed 
maximum value (±1. 7E38) is returned if the dividend is nonzero. 

An "SM" error message is reported if an arithmetic overflow 
occurs. The signed maximum value (±1.7E38) is returned. 

A value of 0 is returned for an overflow. 


N$22, F$ER 
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Purpo se 


To divide a complex number by a real number. 


DAP Calling CALL D$52 

Sequence DAC ARG2 (a real number) 

(Return) 


Method 


This subroutine divides the complex value in the complex 
accumulator (Y) by the real argument, ARG2. 

Y/ARG2 - (A+ B-'I)/ARG2 - A/ARG2+ B : 'T, where Y — A+ B^T 


Data Type of 
Arguments and 
Results 


< implicit complex argument> /< real argument> < complex 
re suit > 


-Other Routines F$AT, H$55, SUB$, L$22, D$22, H$22, L$55 

Used 
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Purpo se 


To divide two complex numbers. 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Re suits 


Other Routines 
U sed 


CALL D$5 5 

DAC ARG2 (complex divisor) 
(Return) 


The following algorithm is used to compute the operation of 
dividing two complex numbers. The contents of the complex 
accumulator (X) are divided by the contents of ARG2 (Y). 

X/Y = (A+ B*I)/(M+ N*I) 

where X - A+ B*I and Y - M+ N'‘T 

= (A+ B*I)*(M-N*I)/ (M+ N*I)*(M-N*I) 

= (A+ B*I)*(M-N*I)/ (M**2+ N**2) 

- (A*M+B*N+ B*M*I-A*N*I)/(M**2+N**2) 

= (A*M+ B*N)/(M**2+ N**2)+ (B*M*I- A*N*I) / (M**2+ N**2) 

= (A*M+ B*N) / (M**2+ N**2)+ (I*(B*M-A*N))/(M**2+ N**2) 

< implicit complex argument> / < complex argument> < complex 
re suit > 


F$AT, H$55, SUB$, L$22, M$22, H$22, A$22, D$22, S$22, 
N$22, L$55 
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Purpo se 

DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To divide a double-precision number by a real number. 


CALL D$62 

DAC ARG2 (a double -precision number) 
(Return) 


This subroutine calls DBLE to convert the real divisor (ARG2 ) 
to a double -precision number and calls the double-precision 
divide routine (D$66). 


< implicit double -precision argument > / < real argument > — 

< double -precision result> 


F$AT, H$66, DBLE, L$66, D$66 
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Purpose 


To divide normalized double -precision numbers. 


See A$6b. 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To calculate the value of an integer raised to an integer power. 
(This subroutine has a high-speed version, E$11X.) 


CALL E$1 1 

DAC ARG2 (the integer exponent) 
(Return) 


The implicit integer argument in the A-register and the integer 
exponent, ARG2, are first examined for the combinations listed 
below. If one of these combinations is found, the answer is loaded 
in the A-register for return to the calling program. 


Value in A -Register Exponent Answer 

I 0 1 

0 0 1 

0 - +32767 

0 +0 

1 J 1 

-1 even 1 

-1 odd -1 

1 - 0 


Otherwise, the value of the expression is calculated and returned 
in the A-register. The maximum or minimum value computed 
may not exceed +32, 767 or -32, 768. 


< implicit integer argument> ** < integer argument> — ► < integer 
re suit > 


The error message "II" is reported and +32,767 is returned if 
overflow occurs or if I = 0 and J is negative (1/0). The value 
-32,768 is returned if 1^-2, J is odd, and overflow occurs. 


ARG$, M$ll, F$ER 


5-2 9 


AM74 




Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 

Other Routines 
Used 


To calculate the value of an integer raised to an integer power. 
(This subroutine requires the High-Speed Arithmetic Option.) 


CALL E$ 1 1 X (or E$ll) 

DAC J (the integer exponent) 

(Return) 


See "Method" for E$ll. 


< implicit integer argument> **< integer argument> — < integer 
re suit > 


See "Error Messages" for E$ll. 


ARG$, F$ER 
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Purpose 


To calculate the value of a real number raised to an integer power. 


DAP Calling CALL E$21 

Sequence DAC ARG2 (the integer exponent) 

(Return) 


Method A**ARG2 is evaluated by multiplying A by itself ARG2-1 times. 

The sign is determined by the sign of the number in the A- and 
B- registers and whether I is odd or even. 

Data Type of . < implicit real argument> **< integer argument> — ► < real result> 

Arguments and 

Results 


Other Routines ARG$, M$22, D$22 

Used 
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Purpose 

DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Re suit s 


Other Routines 
Used 


To calculate the value of a real argument raised to a real power. 


CALL E$22 

DAC ARC 2 (the real exponent) 

( R eturn) 

X**ARG2 is evaluated as e**(ARG2*log(X)). 

< implicit real argument> ** < real argument > — - < real result> 


ARGS, ALOG, M$2 2 , EXP 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Re suits 


Other Routines 
Used 


To calculate the value of a real number raised to a double -precision 
power. 


CALL E$26 

DAC ARG2 (the double-precision exponent) 
(Return) 


B#*ARG2 is evaluated as e**ARG2*log(B)). 


< implicit real argument> ** < double -precision argument> — 

< double-precision result> 


F$AT, C$26, H$6b, DLOG, M$66, DEXP 
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Purpose 


To calculate the value of a complex quantity raised to an integer 
power. 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL E$51 

DAC ARG1 (the integer exponent) 
(Return) 


The number in the complex accumulator is multiplied by itself 
ARG1-1 times. 


<implicit complex argument> ** < integer argument> < complex 
result > 


F$AT, H$55, IABS, L$55, M$55, D$55 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To calculate the value of a double-precision number raised to an 
integer power. 

CALL E$6 1 

DAC ARG2 (the integer exponent) 

(Return) 


This routine checks for an even-numbered exponent, squares the 
number in the double-precision accumulator, and divides the integer 
argument (the exponent) by 2 until the exponent divided by 2 = 1. 

If the exponent is odd, the computed value (D^'l) is multiplied by 
the original double-precision number before exiting. 


< implicit double -precision argument> ** < integer argument> — *- 

< double -precision result> 


F$AT, H$66, L$66, D$66, D$ll, M$ll, M$66 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 

Other Routines 
Used 


To calculate the value of the number in the double -precision accumu- 
lator raised to a real power. 


CALL E$62 

DAC ARG2 (the real exponent) 
(Return) 


B**ARG2 is evaluated as e**(ARG2*DLOG(B)), where B = the 
contents of the double-precision accumulator. 


< implicit double -precision argument> ** < real argument> — ► 

< double -precision result> 


F$AT, H$66, DLOG, M$62, DEXP 
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Purpose 


To calculate the value of a double-precision value raised to 
double-precision result. 


DAP Calling CALL E$66 

-?.ffl^g. nce DAC ARG2 (the double-precision exponent) 

(Return) 


Method 


B**ARG2 is evaluated as e**(ARG2*LOG(B)), where B = the 

contents of the double-precision accumulator. 


Data Type of < implicit double -precision argument > #< double -precision 

Arguments and < double -precision result> 

Results 


Other Routines F$AT, H$66, DLOG, M$66, DEXP 

Used 


argument 
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Purpose 

DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


To store (hold) the contents of the A- and B -registers in memory. 
CALL H$22 

DAC ARG1 (location in which the contents of the A- and B- 

(Return) registers are to be stored) 


The contents of memory at the location specified by the argument 
address, ARG1 , are replaced by the contents of the A- and B- 
registers. The contents of the A- and B -registers remain un- 
changed. 


This subroutine stores a real number in the argument address. 


ARG$ 
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Purpose 


To hold (store) the contents of the complex accumulator in memory 


DAP Calling 
Sequence 


CALL H$55 

DAC ARG1 (location in which the contents of the complex 
(Return) accumulator are to be stored) 


Method The contents of memory at the location specified by the argument 

address, ARG1, are replaced by the contents of the complex accumu- 
lator. The contents of the accumulator remain unchanged. 


Data Type of This subroutine stores a complex number in the argument address 

Arguments and 

Results 


Other Routines ARG$, AC1, AC2, AC3, AC4 

Used 
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Purpose 


To hold (store) the contents of the double -precision accumulator 
in memory. 


DAP Calling 
Sequence 


CALL H$66 

DAC ARG1 (location in which the contents of the double - 

(Return) precision accumulator are to be stored) 


Method The contents of memory specified by the argument address, ARG1 

are replaced by the contents of the double-precision accumulator. 
The contents of the accumulator are unchanged. 


Data Types of 
Arguments and 
Results 


This subroutine stores a double -precision number in the argument 
address . 


Other Routines ARG$, AC1, AC2, AC3 

U sed 
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Purpose 

DAP Calling 
Sequence 

Method 

Other Routines 
Used 


To load a real number into the A- and B- registers. 


CALL L$22 or CALL REAL 
DAC ARG1 (a real number) 

(Return) 


This subroutine calls ARG$ to place the address of the argument, 
ARG1, into the index register. ARG1 is then loaded into the A- 
and B- registers. 


ARG$ 
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Purpose 


To form an INCLUSIVE OR from memory with the value in the 
A- register. 


DAP Calling 
Sequence 


CALL L$33 

DAC ARG1 (an integer value) 
(Return) 


Method 


The value in the A- register is EXCLUSIVELY ORed, ANDed, 
EXCLUSIVELY ORed again with the argument, ARG1 . 





CALL L$55 

DAC ARG1 (a complex number) 
(Return) 


Method 


This subroutine calls ARG$ to place the address of the argument, 
ARG1, into the index register. ARG1 is then loaded into the complex 
accumulator. 


Other Routines ARG$, AC 1 , AC2, AC 3, AC4 

Used 
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Purpose 


To load a double-precision number into the double-precision 
accumulator. 


DAP Calling 
Sequence 


CALL L$66 

DAC ARG1 (a double-precision number) 
(Return) 


Method This subroutine calls ARG$ to place the address of the argument, 

ARG1 , into the index register. ARG1 is then loaded into the 
double-precision accumulator. 


Other Routines ARG$, AC 1 , AC2, AC 3 

Used 
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M$11 


Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To multiply two integers. (This subroutine has a high-speed version, 
M$1 IX. ) 

CALL M$1 1 

DAC ARG2 (integer multiplier) 

(Return) 


This subroutine multiplies the value in the A- register by the integer 
argument, ARG2. If either or both are negative, a sign counter is 
incremented and the negative value(s) are made positive. The multi- 
plier, ARG2, is loaded into the B- register and shifted to place the 
low- order bit of the multiplier in the C- register. The C-bit is 
tested and if it is set, the multiplicand is added to the A- register. 
The A- and B- registers are shifted together 1 bit, with the new low- 
order bit going into the C- register, and so forth, for 16 shifts. 

When these right shifts are completed, the bits are shifted back into 
the A- register, one at a time, checking for overflow. The positive 
or negative result is returned in the A- register. 


< implicit integer argument> *< integer argument? -► < integer result> 


When an over /underflow occurs, the error message "IM" is reported. 
The subroutine returns with +32,767 in the A-register if the answer 
is positive, or -32, 768 if it is negative. 


ARG$, F$ER 
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M$11X 


Purpose 


To multiply two integers, (This subroutine requires the High-Speed 
Arithmetic Option, ) 


DAP Calling 
Sequence 


CALL M$ 1 1 X (or M$1 1 ) 

DAC ARG2 (an integer value) 
(Return) 


Method This subroutine multiplies the value in the A- register by ARG2. 

" “ The result is then examined for over/ underflow (see "Error 

Messages"). If the result is in the proper range, the signed result 
is returned to the calling program in the A- register. 


Data Type of <implicit integer argument > * <integer argument > — cinteger result > 

Arguments and 

Results 


Error Messages See "Error Messages" for M$ll. 


Other Routines ARG$, F$ER 

Used 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Error Messages 


Other Routines 
Used 


To multiply or divide two real numbers. (This subroutine has a 
high-speed version, M$22X. ) 


CALL M$22 
DAC ARG2 
(Return) 


(or D$22) 
(multiplier 
or divisor) 


The multiplicand (M$22) or dividend 
(D$22) must be in the A- and 13- 
registers. The sign, exponent, and 
most significant bits will be in the 
B- register. 


X*Y = (X*2**B)*(Y*2**C), where X = the value in the A- and B- 

register s 

Y = ARG2 

= ABS(X)*ABS(Y)*2**(B+C) 

ABS(X)*ABS(Y) = X(1)*Y(1)+(X(1)*Y(2)+X(2)*Y(1))*2**- 15 
The most significant part of the product is H(X(1)*Y(1)) and 
the least significant part is L(X(1)*Y(1))+H(H(1)*Y(2))+H(X(2)* 
Y(l))*2**-.15. 

Newton's method for 1/Y is R{1) = R(0)*(2- R(0)* Y), where 
R(0) = 1/H(Y), H(Y) being the high-order 15 bits of Y. 

X(l/Y) = X*R(1) = X*R(0)*(2-R(0)*Y). 


< implicit real argument> *< real argument> — ► < real result> 


Multiplication - If there is underflow, a value of zero is returned 
with no error message. 

If there is overflow, an "SM" error message is reported and the 
maximum value ((2**128)- 1) is returned in the A- and B-registers. 

Division - If division by zero is attempted, a "DZ" error message 
is reported and the result in the A- and B-registers is undefined. 

If the divisor is unnormalized, an "SD" error message is reported 
and the result in the A- and B-registers is undefined. 


N$22, ARG$, F$ER 
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Purpose 


To multiply two real numbers. 


DAP Calling 
Sequence 


Method 


Dat a Type of 
Arguments and 
Re suits 


CALL M$22 

DAC ARG2 (a real number) 

(Return) 

X* Y = ( X* 2 * * B) * ( Y * 2 * * C) , where X = the value in the A- and 13- 

registers 

Y = ARG2 

= ABS(X)*ABS(Y)*2**(BC) 

ABS(X) : '-ABS(Y) = X( 1)* Y( 1)*(X( 1)* Y(2)+X(2)* Y( 1))*2** - 1 5 
The most significant part of the product is H((X( 1)*Y{ 1)) and the 
least significant part is L(X(1)*Y(1))+H(H(1) : '-'Y(2))+H(X(2)''-Y(1)) 

15, where H(X(1)-Y(1)) is the most significant part of the 
product X ( 1 ) * Y ( 1 ) and L(X(1)*Y(1)) is the least significant part 
of that product . 

< implicit real argument> * < real argument> -► < real result> 


Error Messages Un derflow - A value of zero is returned with no error message. 

— — - i ~~~ Overflow - An "SM" error message is reported and a signed maxi- 

mum value (±1.7E38) is returned. 


Other Routines F$ER 

Used 
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Purpose 


To multiply a complex number by a real number. 


DAP Calling CALL M$52 

Sequence DAC ARG2 (a real number) 

(Return) 


Method 


Y*X = (A+B*I)*X = A*X+(B*X)*I 
where Y = A+B*I (in the complex accumulator) 
X = ARG2 


Data Type of < implicit complex argument> * < real argument> — ► < complex result> 

Arguments and 

Results 


Other Routines F$AT, H$55, SUB$, L$22, M$22, H$22, L$55 

Used 
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Purpose 


To multiply complex numbers. 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 

Used 


CALL M$55 

DAC ARG2 (a complex value) 
(Return) 


This routine multiplies the contents of the complex accumulator (X) 
by the value in ARG2 (Y). 

X* Y = (A+B*I) (M+N*I) = A*M-B*N+(A*N+B*M)*I 
where X = A+B*I and Y = M+N*I. 


< implicit complex argument> * < complex argument> — ► < complex 
re suit > 


F$AT, H$55, SUB$, L$22, M$22, H$22, S$22, N$22, A$22, L$55 
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Purpose 


To multiply a double-precision number by a real number. 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL M$62 

DAC ARG2 (real multiplier) 
(Return) 


This subroutine calls DBLE to convert the real multiplier to a 
double-precision number and calls the double-precision multiply 
routine (M$66). 


< implicit double -precision argument? * < real argument? — *- 

< double -precision result> 


F$AT, H$66, DBLE, M$66 
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Purpose 


o multiply normalized, double-precision numbers. 


See A$66. 
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Purpose 


To determine the TWOs complement of a real number. 


DAP Calling 
Sequence 


CALL N$22 

DAC ARG1 (a real number) 
(Return) 


Method 


The C-bit is preset on entrance to this routine to provide a true TWOs 
complement if the low- order word is found to be zero. The C-bit is 
reset when this is not the case, and the A- and B- registers are TWOs 
complemented normally. 


Data Type of The TWOs complement of the real argument is computed and the 

- gfg uments ^d routine exits with the the real result in the A- and B- registers. 

Results 6 
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To obtain the complement of a logical value. 



CALL N$33 
(Return) 

The least significant bit of the argument in the A- register is 
logically complemented, changing its value from true to false 
(1 to 0) or false to true (0 to 1). 
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To negate a complex quantity. 


Purpose 


DAP Calling 
Sequence 


CALL N$55 
(Return) 


Method 


The signs of the real part and the complex part of the complex 
number are negated. The result is in the complex accumulator. 


Data Type of 
Arguments and 
Results 


The complex argument is negated and the subroutine exits, with the 
complex result in the complex accumulator. 


Other Routines 
Used 


H$55, SUB$, L$22, N$22, H$22, L$55 
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Purpose 


To negate a double-precision number. 


DAP Calling 
Sequence 


CALL N$66 
(Return) 


Method 


This subroutine negates the value in the double-precision accumu- 
lator. The double-precision word is effectively TWOs complemente 
as follows: 

1. The lowest order word, AC3, word 3, is tested for zero. If it 
is not zero, the word is TWOs complemented. If it is zero, the 

C-bit is set. 

2 AC2, word 2, is tested for zero. If it is not zero, the word is 
ONEs complemented and the C-bit is added. If it is zero and the 
C-bit is set, no action is taken. If the C-bit is not set, the word 
is ONEs complemented. 

3. AC1, word 1, is ONEs complemented, and the C-bit, if set, is 

added. The negated result is left in the double-precision accumu- 
lator. 


Data Type of The double-precision argument is negated and^the routine exits with v 

Arguments and a double-precision result in AC 1 , AC2, and AC 3. 

Results 


Other Routines AC1, AC2, AC3 

Used 
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Purpose 


To subtract real numbers. (This subroutine has 
a high-speed version, S$22X. ) 


See A$22 . 
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Purpose 


See A$22X. 


To subtract real numbers. (This subroutine requires the 
High-Speed Arithmetic Option. ) 
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Purpose 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


S$52 

To subtract a real number from a complex number to obtain a com- 
plex result. 


CALL S$52 

DAC ARG2 (a real number) 

(Return) 

Y-X = A+B*I-X = (A-X)+B*I 
where Y = A+B*I, X = ARG2 


< implicit complex argument> -<real argument? -► < complex result? 


F$AT, H$55, L$22, S$22, H$22, L$55 
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Purpose 


To subtract two complex numbers. 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
R e suits 


Other Routines 
U sed 


CALL S$55 

DAC ARG2 (the complex subtrahend) 

(Return) 


This subroutine subtracts ARG2( Y) from the value in the complex 
accumulator (X): 

X - Y = (A+B*I) - (M+N*I) = A*M+(B*N)*I 
where X = A+B*I. Y = M+MI 


< implicit complex argument> — <complex argument> < complex 

re suit > 


F$AT, H$55, SUB$, L$22, S$22, N$22, H$22, L$55 
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Purpose 


To subtract a real argument from a double -precision number. 


DAP Calling 
Sequence 


Method 


Data Type of 
Arguments and 
Results 


Other Routines 
Used 


CALL S$62 

DAC ARG2 (a real number) 

(Return) 


This subroutine calls DBLE to convert the real argument to a double - 
precision number and enters the double -precision subtraction 
routine (S$66). 


< implicit double -precision argument> — < real argument> 

< double -precision result> 


F$AT, H$66, DBLE, S$66, N$66 
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Purpose 


To subtract normalized, double-precision numbers. 


See A$66. 
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Purpose 


To convert a double-precision number to a real number. 


See C$62 . 
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Purpose 


To calculate the address of a referenced array element or to calcu- 
late the array size. 


DAP Calling 
Sequence 


CALL SUB$ 

DAC or DAC* ARRAY TABLE1 
DAC or DAC* SUBSCRIPT 1 
DAC or DAC* SUBSCRIPT 2 


DAC or DAC* SUBSCRIPT N 
(Return) 


or 

CALL SIZ$ 

DAC or DAC* ARRAY TABLE2 
(Return) 


Method 


ARRAY TABLE1 Layout 


DAC or DAC* ARRAY 

OCT L (number of words per array element) 

DEC DIMENSION 1 

DEC DIMENSION 2 

DEC DIMENSION N 

OCT O (end of dimension list) 

ARRAY TABLE2 Layout 

DAC or DAC* ARRAY 
OCT KEY 

OCT or DAC* DIMENSION 1 
OCT or DAC* DIMENSION 2 


OCT or DAC* DIMENSION N or OCT ARRAY SIZE 
or omitted 

The KEY bit pattern is CVDDDDDDDDDDDLLL, where 

C = 0 - no array bounds checking 
C = 1 - array bounds checking 

V = 0 - last word of array table is array size 

V = 1 - last word of array table is dimension 

If C = 0 and V = 0, the last dimension word of the array table is 
omitted. 

D = dimensionality - limited to 2047 
L = number of words per array element 
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Note that L is determined by the data type of the array as follows: 


Data Type of Array 

L = 1 - integer or logical 

2 - real 

3 - double-precision 

4 - complex 

Let S denote the array starting address, L the number of words per 
array element, S(I) the Ith subscript value, and D(I) the Ith dimension 
for an N-dimensional array A where N > 1. 

The address of the array element A(S(1), S(2), , . . S(N)) is given by 
S+ L{..,, (S(N)-1)*D(N-1) +...+ (S(2)-1)*D(1) + (S(l)-l)) 


Error Messages The error message "AO" (array overflow) is reported if the array 

element referenced is outside the bounds of the array. Only the final 
array element referenced is checked for legality, not individual sub- 
script values. 

Other Routines M$ll, F$ER 

Used 
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Purpose 


To clear (zero- out) the exponent of the variable in the double- 
precision accumulator. 


DAP Calling 
Sequence 


CALL Z$80 
(Return) 


Method 


Extract the value in AC1 and replace the characteristic (base 2) in 
bits 2-9 with zeros. 


Other Routines AC 1 

Used 
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SECTION VI 


RUN-TIME AND CONTROL SUBROUTINES 

This section describes the routines which: control input and output by selecting and 
activating the proper device drivers; provide buffers; and edit and trace all I/O. 
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Purpose 

To transfer 

an array 

from or to the input or output data list, 

DAP Calling 

CALL 

F$AR 

For use with DAP 

Sequence 

OCT 

m 

Number of words in array 

DAC 

a 

Location of first word 


(Return) 




or 


• 


CALL 

F$Lx 

As used by the compiler 


OCT 

m 

x = 1 for integer 


DAC 

a 

= 2 for real 


(Return) 


= 3 for logical 


= 5 for complex 
= 6 for double precision 


FORTRAN 

Reference 


DIMENSION 1(5) 

T>TT* ATI Kof of 


WRITE (x, f) list 


l f — TT 1 n T? r T’ T7 A TVT pfofpmpnf 


Method Whenever the data list requires data from the internal source 

or data to be stored in the internal source, exit is made from 
F$10 to the next location in the data list. Elements of the data 
list may be variables, subscripted variables, or array names. 
Mode may be integer, real, logical, complex, or double 
precision. For easier data transmission, all list elements 
are assumed to be arrays; the mode is determined by 
appropriate format descriptors. For each item in the list, 
the three -word calling sequence (above) is generated. 


Other Routines F$10, F$CB, F$ER 

Used 
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Purpose 


DAP Calling 
Sequence 


FORTRAN 

References 


Method 


Data Type 
of Arguments 


To transfer a variable number of arguments from the 
calling routine to the called routine. 


DAC 

CALL 

DEC 

DAC 


** Entry point 

F$AT 

n Number of arguments to be transferred 

ARG1 Address at which first argument is stored 


DAC ARGn 

(Return) 


Address at which last argument is stored 


CALL SUB 1 (ARG1 , ARG2, ARGn) 

Where SUB1 is any subroutine and ARG1 
through ARGn are any constants, variables, 
arrays, etc. 


When arguments are to be transferred from a calling routine 
to a subroutine, a call to F$AT is generated by the compiler. 
The number of arguments specified by the first pseudo - 
operation following the call are transferred from the calling 
routine to the subroutine. All levels of indirect addressing 
are removed before an argument is transferred. ARG1 is 
the beginning location of the block into which the arguments 
are to be placed. 


Arguments are direct relative addresses. 
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Purpose 


To connect the calling program with the magnetic tape 
rewind routine. 


PAP Calling 
Sequence 


CALL F$Bx x = 5, 6 , 7, 8 , 9 or the previously defined 

variable n (n = 5, 6 , 7, 8 , or 9) 


FORTRAN REWIND x x = 5 , 6 , 7, 8 , 9 or the variable n 

Reference 


Method This routine converts the logical magnetic tape unit number 

— £, 7 ( g or 9 ) to the corresponding physical magnetic tape 

unit number (1,2, 3, 4, or 5) and then calls the REWIND 
routine to rewind the tape on that unit to the beginning of tape. 


Uther Routines C $MR 

Used 
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Purpose 

DAP Calling 
Sequence 

FORTRAN 

Reference 


Method 


Other Routines 
Used 


To close the buffers used for input or output. 


CALL 

F$CB 

READ 

(n, m) list 

READ 

(n, m) 

READ 

(n) list 

READ 

(n) 

WRITE 

(n, m) list 

WRITE 

(n, m) 

WRITE 

(n) list 


n = device number; m = format 
statement number 


At the end of the data list from a READ or WRITE statement, 
or when the format statement is exhausted (non-list), a call 
is issued to F$CB to close the buffer. The address of the 
buffer was determined by F$IO. F$CB checks for I/O mode 
and immediately closes the buffer if mode is input. If the 
buffer is formatted output, F$CB fills the remainder of the 
buffer with up to 134 spaces. If output is binary, the end of 
the buffer is filled with up to 120 zeros. 


F$IO 
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Purpose 


To control the writing of an end-of-file mark on magnetic 
tape and a STOP code as an end-of-file on paper tape punch. 


DAP Calling 
Sequence 


CALL F$D2 

CALL F$Dx 

CALL F$Dn 


2 = paper tape punch 

x = logical tape unit number 5 through 9 
n = dummy device number and the A- 
register contains the value 2, 5, 6, 

7, 8, or 9 


FORTRAN END FILE x x = 2, 5, 6, 7, 8, 9 or the variable n 

References 

Method This routine converts the logical magnetic tape unit number 

(5, 6, 7, 8, or 9) to the corresponding physical magnetic tape 
unit number (1, 2, 3, 4, or 5). It then calls the driver to 
write an end -of-file mark on the specified magnetic tape. 

If x is 2, it calls the driver to punch a STOP code on the 
high-speed paper tape punch. 


Other Routines G$ME, 0$PS 

Used 
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Purpose 


DAP Calling 
Sequence 


FORTRAN 

References 

Method 


Data Type 
of Argument 

Other Routines 
Used 


To cause a mnemonic error indicator to be typed on the ASR-33 
when an object-time error is encountered in a specified routine. 


CALL F$ER ARG1 is the address of the indicator to 

DAC ARG1 be typed. The routine types the error 

indicator and halts if Sense Switch 3 
is not set. Pressing START after the 
halt causes the program to continue. 

If Sense Switch 3 is set, F$ER exits 
with no typeout and no halt. 


CALL F$ER (2Hxx) xx = two ASCII characters to be typed 


The location of the object-time error mnemonic indicator 
is extracted from its relocatable address, AC5. AC 5 is 
then loaded with the error mnemonic indicator and Sense 
Switch 3 is tested. If the switch is not set, a carriage 
return and line feed are issued. The error mnemonic is 
then printed, and the routine halts. If START is pressed 
at this point, a normal return is made to the calling routine. 
If Sense Switch 3 is set, return is made immediately to 
calling routine. 


The argument is the address of any two ASCII characters. 


AC 5, F$HT 
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Pu r po s e 


To control back spacing of a record on magnetic tape. 


DAP Calling CALL F$Fx x = 5, 6, 7, 8, 9 or the previously defined 

Sequence variable n. 


FORTRAN BACKSPACE x x = 5, 6, 7, 8, 9 or the variable n. 

References 


Method This routine converts the logical magnetic tape unit number 

(5, 6, 7, 8, or 9) to the corresponding physical magnetic tape 
unit number (1,2, 3,4, or 5) and then calls the driver to back 
space one record on the specified magnetic tape. 


Error 

Message 


An error message BF is reported if an end-of-file is 

encountered. 


Other Routines C$BR, F$ER 

Used 
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Purpose 

To process FORTRAN run-time assigned GO TO statements. 

DAP Calling 
Sequence 

LDA PTR 

CALL F$GA 

DEC n 

DAC SI 

DAC S2 

Transfer address in A-register 

Number of statements in list 

Address of first statement 

Address of second statement 


DAC Sn 

Address of last statement 

FORTRAN 

Reference 

ASSIGN J TO I 

GO TO I, (Kl, K2, . . . 

J = statement number 
i Kn) I = integer variable name 

Ks = statement numbers 

Method 

This routine checks the address passed in the A-register against 
the statement address list that follows the call. If the address 
is found in that list, control passes to that statement. If not, 
a GO error is reported. No recovery from this error is 
possible. If the statement address list is empty (n = 0), the 
address is not checked. 

Data Type 
of Argument 

An address is passed to this routine in the A-register. Con- 
trol is passed to the statement number at that address. 

Other Routines 

Used 

F$ER 
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Purpose 


To process FORTRAN run-time computed GO TO statements. 


DAP Calling 
Sequence 


LDA 

PTR 

Index to 

statement list 

CALL 

F$GC 



DEC 

n 

Number 

of statements in list 

DAC 

SI 

Address 

of first statement 

DAC 

S2 

Address 

of second statement 


DAC Sn 


Address of last statement in list 


FORTRAN 

Reference 


GO TO (Kl, K2 , . . . Kn), I I = integer value in the range 1 to n 

Ks = statement numbers 


Method The integer variable I or the content of PTR is treated as an 

index number; F$GC uses it to select the statement number 
from that position in the calling sequence. For example, 
if I = 3 then control is shifted to statement K3 or S3 in the 
above examples. If the index (I) is < 1 or >n, the computed 
GO TO statement is treated as a NOP. 


Data Type 
of Argument 


An integer value is passed to this routine in the A-register. 
Control is passed to the statement at the computed address. 
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F$HT 


Purpose To cause the computer to stop and print PA if a PAUSE 

statement has been encountered, or to print ST if a STOP 
statement has been encountered. 


DAP Calling 
Sequence 

CALL 

DAC 

CALL 

DAC 

F$HT 

'151724 Octal notation for ST 

F$HT 

'150301 Octal notation for PA 

FORTRAN 

Reference 

STOP or 

PAUSE 


Method 

This calling sequence is generated by the compiler when the 
STOP or PAUSE verb is encountered. The mnemonic ST or 

PA is placed in relocatable address AC 5 and printed on the 

ASR. The A-register is restored and the program halts. 

Return to the calling program may be made by pressing START 

Data Type 
of Argument 

The binary < 

equivalent to the specified ASCII characters. 

Other Routines 

Used 

AC5 
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FSIO 


Purpose 


DAP Calling 
Sequence 


FORTRAN 

References 


Method 


Other Routines 
Used 


To perform input/output conversion, to edit input/output 
information, to accommodate the appropriate input /output 
device, and to provide buffer s . 


DAC 

CALL 

DAC 


a a = location of the format list from the 

F$IO READ statement 

BUF BUF = buffer location 


DAC* a The calling sequence for a FORTRAN 

CALL F$IO WRITE with a and BUF as above 

DAC BUF 


READ (n, x) list n = device number 

WRITE (n, x) list x = format statement number 

x FORMAT ( ) 


A FORTRAN READ / WRIT E statement starts with a device 
number and a reference to a format statement, followed 
by an optional argument list. The first instruction generated 
by the READ/ WRITE statement is a coupling to the appro- 
priate device driver. The device driver then calls on 
F$IO passing the location of the format list, setting the entry 
location for the device driver and setting a flag to indicate 
input or output. F$IO then interprets the format list, 
character -by-character, taking whatever actions are 
required. Whenever data is required from or is to be stored 
in the internal source, exit is made from F$IO to the next 
location in the data list. 


F$AR, F$CB, F$ER 
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F$R1 


Purpose 

DAP Calling 
Sequence 

FORTRAN 

References 

Method 


Data Type 

Other Routines 
Used 


To control the typewriter keyboard input routine. 


CALL F$R1 

DA C n n = location of the format descriptor list 

(Return) 


READ(1, f) list f - FORTRAN statement number 

READ(1, f) 


This subroutine connects the calling program with the 
I/O control subroutine (F$IO). Included in F$R1 is the driving 
logic needed to transmit input from the typewriter keyboard. 
When F$IO is called, the location of the format descriptor 
list (if any), the entry location of the driver subroutine, and 
a flag indicating input are transmitted. 

Whenever the F$IO subroutine requires data, return is made 
to the driver input entry of this subroutine, at which time up 
to 120 characters (terminated by a carriage return) are 
entered into the input buffer. 


Information is in ASCII format. 


F$IO, I$AA 
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Purpose 


To control the paper tape input routine. 


DAP Calling 
Sequence 


CALL F$R2 

DAC n 

(Return) 


Location of the format descriptor list 
(or 00000 if input is in binary format) 


FORTRAN 

References 


READ (2, f) list 
READ (2, f) 
READ (2) list 
READ (2) 


Formatted READ where f is a format 
statement number 
Unformatted paper tape read 


Method This subroutine connects the calling program with the I/O 

control subroutine (F$IO). Included in F$R2 is the driving 
logic needed to transmit input from the paper tape reader. 

When F$IO is called the location of the format descriptor 

(; .c ~ v— 1r>c>of^n of f V. ^ A v'i r- cnhrniitinp. and 

xxoc > y > j 

a flag indicating input are transmitted. 

Whenever the F$IO subroutine requires data, it calls on the 
F$R2 driver subroutine, which assembles data into a 60-word 
buffer, three characters per word, if the input mode is 
binary. If the input mode is BCD, 80 characters are 
assembled, two per word, into a 40 -word buffer. A carriage 
return character is replaced with as many blanks as needed 
to fill the rest of the input buffer. A tab character is 
replaced with as many blanks as needed to reach the next 
predetermined tab position. 

A STOP code read in either binary or BCD mode causes the 
characters ST to be typed followed by a halt. Press START 
to continue. 


Data Type 


Information is in ASCII if formatted, or in binary if 
unformatted. 


Other Routines F$IO, I$PA, I$PB 

Used 
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Purpose 

DAP Calling 
Sequence 

FORTRAN 

References 

Method 


Data Type 

Other Routines 


To control the card input routine. 


CALL F$R3 

DAC n 

(Return) 


Location of format descriptor list 
(00000 if input is binary) 


READ (3, f) list f - FORTRAN Statement number 

READ (3, f) 

READ(3) list Unformatted read from a punch card 

READ (3) 


This subroutine connects the calling program to the I/O 
control subroutine F$IO. Included in this subroutine is the 
driving logic required to input from the card reader. When 
F$IO is called, this subroutine transmits the location of the 
driver subroutine and a flag indicating input. 

When the F$IO subroutine requires data, it calls upon the 
F$R3 driver subroutine, which assembles 80 characters 
(two per word) into a 40 -word buffer, if mode is Hollerith 
(formatted), or into a 60 -word buffer in column binary 
format (three words for every four columns) if binary. 
Return is made to F$IO where the buffer is processed. 


Information is in Hollerith if formatted, or in column 
binary if unformatted. 


F$IO, I$CA, I$CB 


Used 



F$R5-9 

Purpose 

DAP Calling 
Sequence 

FORTRAN 

References 

Method 


Data Type 


Other Routines 
Used 


To control reading of magnetic tape. 


CALL F$Rx 

DAC n 

(Return) 


x = 5 , 6, 7, 8, or 9 

Location of the format descriptor list 
if formatted, or zero if unformatted 


READ(x,f) list 
READ (x, f) 
READ (x) list 
READ(x) 


f = FORTRAN statement number and 
x=5,6,7,8, or 9 

Unformatted read where x = 5,6, 7, 8, or 9 


This subroutine connects the calling program with I/O control 
routine (F$IO) and standard magnetic tape routines. 

When F$IO is called, the format descriptor list and a flag 
indicating input are transmitted. 

When the F$IO routine needs a buffer of data, it calls this 
driver, which in turn calls the appropriate magnetic tape unit 
and conversion routines (for formatted READ). The number 
of words read is 60, equivalent to 120 characters. 

The appropriate magnetic tape units are physical magnetic tape 
units 1 through 5 corresponding respectively to logical tape 
units, numbers 5 through 9. 


Information is in ASCII if formatted, or in binary if 
unformatted. 


F$IO, I$MA, I$MC , C$6T08 


6-16 


AM74 



Purpose 

DAP Calling 
Sequence 

FORTRAN 

References 


Method 


Data Type 


Other Routines 
Used 


To control the input drivers for variable input device numbers. 


LDA d 

CALL F$Rn 

DAC n 

(Return) 


Location of device number 

Location of Format Descriptor List = 00000 
if format is binary 


READ (x, f) list 
READ (x, f ) 
READ(x) list 
READ(x) 


f - FORTRAN statement number; 
x = variable device number 1 through 9 

unformatted read where x = 1, 2, 3, 5, 6, 7 , 8, 9 


The value of d is checked for correct limits and is then used 
to determine the entry position of a Jump Table. The Jump 
Table transfers to the proper F$R subroutine. (Note that 
the entire F$R subroutine must be called into memory along 
with this subroutine, because there is no way of knowing 
in advance which drivers are required. ) 

If d does not equal a number from 1 to 9, the computer 
halts with a 1 in the A-register. The A-register may be changed 
manually to another device number; otherwise, the typewriter 
will be selected as the input device when START is pressed 
to continue processing. 

Other errors, such as parity, end of tape, etc. , cause the 
actions described in the appropriate F$R subroutine. 


Information is in ASCII if formatted, or in binary if 
unformatted. 


F$R1 , F$R2, F$R3, F$R5-9 
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Purpose 




To aid in debugging programs by printing the following: 

• Values of variables or array elements as they are 
being stored 

• Values of the DO parameters as they vary 

• Locations of statement numbers as they are encountered 

• Content of IF statements as they are evaluated. 


DAP Callin g 
Sequenc e 


CALL F$TR 

OCT 

OCT 

OCT 

(Return) 


The first three bits of the first argument 
are : 

000 = 0 = Statement number 

001 = 1 = Integer 
010 = 2 = Real 
011=3= Logical 

100 = 4 (Not Used) 

101 = 9 r Como] ex 

110 = 6 = Double precision 
The mode of an IF statement depends 
on the mode of the expression being 
evaluated. The remaining 13 bits of 
ARG1 and all the ARG2 and ARG3 are 
the name of the variable, the statement 
number, or blank. 


FORTRAN 

References 


TRACE xl,x2, . . . xn Each x represents a variable or array 

name ; 

TRACE n n represents a statement number 


Method 


The A-register is stored in location zero (index register). 

The location of the argument is extracted and is interchanged 
with the index register, restoring the A-register. Sense 
Switch 4 is interrogated and if set, return is made to the 
calling program. If Sense Switch 4 is reset, trace mode is 
entered. The A- and B -registers are saved, and the buffer 
pointer is reset. The mode jump address is then set up. 

The 1- to 6-character name is then put in the format state- 
ment buffer. The return address for exit is calculated and 
stored. The 1 - to 6 -character name is interrogated to determine 
whether a statement number or a variable is being processed. 

If a variable name is being processed, an equal sign is stored 
in the format buffer following the name of the variable. An 
indirect jump is then taken, contingent upon whether the mode 
is real, integer, logical, complex, o r double -precision. The 
following formats are moved to the format buffer, depending 
on the mode, to provide maximum representation of the 
variable. 


REAL 

INTEGER 

LOGICAL 

COMPLEX 

DOUBLE PRECISION 


G1 9. 7 

17 

L2 

El 5. 7 ,H, El 5. 7 
D 1 9. 11 
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F$TR cont 


The output device and format location are provided to the 
appropriate I/O driver and communicated to F$IO. Then the 
variable name and equal sign, and its appropriate value, are 
printed. When an IF statement evaluation is being written 
the four characters IF ( ) are printed, followed by an equal 
sign and the value of the expression in the parentheses. 

If a statement number is being printed, the number is 
bracketed by opening and closing parentheses, and is printed. 
Upon completion of each printout, the buffer is closed and 
the A- and B -registers restored. Return is made to the 
calling program. 


Data Type See calling sequence. 

of Arguments 


Other Routines F$W1, F$AR, F$CB, AC1, AC2, AC 3 

Used 
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F$W1 


Pur po s e 

DAP Calling 
Sequence 


fortran 

References 


Method 


Data Type 

Other Routines 
Used 


To control the typewriter output routine. 


CALL F$W1 

DAC n Location of the format descriptor list 

(Return) 

WRITE( 1 , f)list f = FORTRAN statement number 

WRITE ( 1 , f) 

This subroutine connects the calling program with the I/O 
control subroutine (F$IO). Included in this subroutine is the 
driving logic needed to produce output on the typewriter. 

When F$IO is called, the location of the format descriptor 
list (if any), the entry location of the driver subroutine, and 
a flag indicating output are transferred. 

After the F$IO subroutine has generated a buffer full of data 
(72 characters), return is made to the driver output entry 
of this subroutine. At that time, the first character of the 
buffer is analyzed for proper line feed control as follows: 

Blank — Type a single carriage return 

followed by characters (2-72) 

0 — Issue two carriage returns, then 

type characters (2-72) 

1 — Skip to the top of the next page, 

then type characters (2-72) 

+ — No line advance, type characters (2-72) 

Others — Same as blank, except characters 
(1 -72 ) are typed 

This subroutine prints 60 lines per page and skips six lines 
to the top of the next page. The operator should start three 
lines down the first page in order to get correct spacing 
between pages. 


Only ASCII information is processed. 


F$IO, 0$AC , 0$AF 
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Purpose 

DAP Calling 
Sequence 

FORTRAN 

References 

Method 


Data Type 

Other Routines 
Used 



To control the paper tape output routine. 


CALL F$W2 

DAC n 

(Return) 


WRITE (2, f) list 
WRITE(2,f) 
WRITE(2) list 


Location of format descriptor list 
(00000 if output mode is binary) 

f = FORTRAN statement number 
Unformatted WRITE to paper tape punch 


This subroutine connects the calling program to the I/O control 
subroutine F$IO. Included in the subroutine is the driving logic 
required to produce output on the paper tape punch. When F$IO 
is called, this subroutine transmits the location of the format 
descriptor list (if any) the entry location of the driver sub- 
routine, and a flag indicating input. 

After the F$IO subroutine has generated a full buffer of data 
60 words at two characters per word, or 40 words at three 
binary characters per word, return is made to the driver 
output entry of the subroutine. At that time, the buffer is 
punched on tape. 


Information is in ASCII if formatted or in binary if unformatted. 


F$IO, 0$PF, Q$PP, 0$PB 
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Purpose 


To control the card punch routines. 


DAP Calling 
Sequence 


FORTRAN 

References 


Method 


Data Type 

Other Routines 
U sed 

Note 


CALL F$W3 

DAC n 

(Return) 


Location of format descriptor list 
(00000 if output mode is binary) 


WRITE(3, f) list f = format statement number 

WRITE(3,f) 

WRITE(3 ) list Unformatted WRITE to the card punch 


This subroutine connects the calling program (FORTRAN 

Object Program) to the I/O control subroutine, F$IO, and to 

the card punch subroutines. When F$IO is called, F$W3 

transmits the location of the address of the format descriptor 

list (if any), including a flag indicating output mode (DAC* for input), 

and a location for r^^ntrance to F*$>W3. 

After the F$IO subroutine has generated a full buffer of data 
(40 words at two BCD characters per word or 60 binary words), 
return is made to F$W3. The appropriate card punch subroutine 
is called, and a card is punched. 


Information is in Hollerith if formatted, or in column binary 
if unformatted. 


F$1C, C$CH, 0$CB 


The Hollerith information can be either 026 or 029 character 
set, depending on the version of Q$CH selected. 
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F$W4 


Purpose 

DAP Calling 
Sequence 

FORTRAN 

References 

Method 


Data Type 

Other Routines 
Used 


To control the line printer output routine. 


CALL F$W4 

DAC n Location of format descriptor list 

(Return) 

If n is zero, which normally indicates binary 
output, the computer halts. Push START to 
print data in BCD format. 


WRITE(4,f) list f = format statement number 

WRITE(4, f) 


This subroutine connects the calling program to the I/O control 
subroutine F$IO. Included in the subroutine is the driving logic 
required to produce output on the line printer. When F$IO 
is called, this routine transmits the location of the format 
descriptor list (if any), the entry location of the driver sub- 
routine, and a flag indicating output mode. 

After the F$IO subroutine has generated a full buffer of data 
(120 characters), return is made to the driver output entry 
of this subroutine. At that time, the first character of the 
buffer is analyzed for proper line feed control as follows: 

Blank — Advance one line and print characters 
2 through 120. 

0 — Advance two lines and print characters 

2 through 120. 

1 — Advance to top of next page and print 

characters 1 through 120. 

+ — Print characters 1 through 120 without 

advancing line position. 

Others — Advance one line and print characters 
1 through 120. 


Information is in ASCII format. 


F$IO, 0$LP, 0$LO 
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F$W5-9 


Pu r po s e 

DAP Calling 
Sequence 

FORTRAN 

References 

Method 


Data Type 


Other Routines 
Used 


To control writing on magnetic tape. 


CALL F$Wx x = 5, 6, 7, 8, or 9 

DAC n Address of format, if any 


WRITE(x, f) list 
WRITE (x, f) 
WRITE (x) list 


f = FORTRAN statement number 
x = 5, 6, 7, 8, or 9 

Unformatted READ, x = 5, 6, 7, 8, or 9 


This program connects the calling program with the I/O control 
routine (F$IO) and the standard magnetic tape routines. WTien 
F$IO is called, the format descriptor list and a flag indicating 
output mode are transferred. 

When the F$IO routine has a buffer of data to write, it calls 
the driver, which in turn calls the appropriate magnetic tape 

' lid ■ Olivers io i rout* res ■ r c r r c vt ett ~ d T ,\ r Pv.IT HP 

60 words are written, equivalent to 120 characters in either 
mode (formatted or unformatted). 

The appropriate magnetic tape units (physical) are numbers 
1 through 5 corresponding respectively to the logical tape 
units numbers 5 through 9 given for x above. 


Information on the tape is in ASCII if formatted, or in 
binary if unformatted. 


F$IO, 0$MC, 0$MA, C$8T06 
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F$Wn 


Purpose 

DAP Calling 
Sequence 

FORTRAN 

References 


Method 


Data Type 

Other Routines 
Used 


To control the output drivers for variable output device numbers. 


LDA d 

CALL F$Wn 

DAC n 

(Return) 


Location of device number 

Location of format descriptor list 
(00000 if format is binary) 


WRITE(x, f) list 
WRITE (x, f) 
WRITE(x) list 


f - FORTRAN statement number 
x = variable device number 1 through 9 
Unformatted WRITE, where x = 1, 2, 3, 4, 5, 
6,7,8, or 9 


The value of d is checked for correct limits and then used to 
determine the entry position of a Jump Table. The Jump 
Table is then transferred to the proper F$W subroutine. 

(Note that all F$W subroutines must be called into memory 
along with this subroutine, because there is no way of knowing 
in advance which drivers are required. 

If d does not equal a number from 1 through 9, the computer 
halts with a 1 in the A-register. The A-register may be 
changed manually to another device number. Otherwise, 
the typewriter will be selected as the output device when 
START is pressed to continue processing. 

Other errors, such as parity, end of tape, etc., cause the 
actions described in the appropriate F$W subroutine. 


Information is in ASCII if formatted, or in binary if unformatted. 
F$W1, F$W2, F$W3, F$W4, F$W5-9 
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APPENDIX A 
TAPE CONTENTS 


MAGNETIC TAPE 70182805000 - LTCSI S 
(LIBRARY SOURCES CODED IN FORTRAN) 


This tape consists of the individual sources of the following programs in the order listed. 
This tape is one of two distributed and contains that portion of the FORTRAN Library that is 
FORTRAN-coded. 


NUMBER 

NAME 

DOC. NO. 

1 

STMEAN 

70181386000 

2 

STGEOM 

70181387000 

3 

STCORR 

70181388000 

4 

STSPER 

70181389000 

5 

STCRMT 

70181390000 

6 

STMEDT 

70181391000 

7 

ST CHI 2 

70181392000 

8 

STBNPB 

70181394000 

9 

STL NRG 

70181395000 

10 

STCHIS 

70181393000 

11 

STPLRG 

70181396000 

12 

STANVI 

70181397000 

13 

STANV2 

70181398000 

14 

STANVR 

70181399000 

15 

STANVL 

70181400000 

16 

STANVG 

70181401000 

17 

STANVB 

70181402000 

18 

STANVY 

70181403000 

19 

DEFOAD 

70181405000 

20 

DEFOMA 

70181404000 

21 

DEFOHA 

70181406000 

22 

DEFORK 

70181407000 

23 

DESOAD 

70181408000 

24 

DESOMA 

70181409000 

25 

DESOHA 

70181410000 

26 

DESORK 

70181411000 

27 

PLYMUL 

70181416000 

28 

PLYDIV 

70181417000 

29 

PLYINT 

70181420000 

30 

PLYIRT 

70181421000 

31 

PLYEVL 

70181422000 

32 

PLYDEF 

70181423000 

33 

NACPLY 

70181429000 

34 

NAAITK 

70181418000 

35 

NALAGR 

70181419000 

36 

NABAIR 

70181424000 

37 

MATTRS 

70181412000 

38 

MATMUL 

70181413000 

39 

MARITH 

70181414000 
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FILE NUMBER 

NAME 

DOC. NO. 

40 

MATEIG 

70181415000 

41 

MAT I NY 

70181427000 

42 

NAREGU 

70181425000 

43 

NAMULL 

70181426000 

44 

GSEID 

70181428000 

45 

SORT 

70181430000 

46 

SORT 2 

70181431000 

47 

CVPOLR 

70181432000 

48 

E$62 

70180053000 

49 

E$61 

70180052000 

50 

E$26 

70182582000 

51 

E$66 

70180054000 

52 

DSQRT 

70182580000 

53 

DCOS 

70180055000 

54 

DSIN 

70182583000 

55 

DEXP 

70182581000 

56 

DLOGIO 

70180051000 

57 

DLOG 

70182579000 

58 

DLOG2 

70182914000 

59 

DATAN2 

70180056000 

60 

DATAN 

70182584000 

61 

DMOD 

/0182b»»UUU 

62 

DSIGN 

70182589000 

63 

DABS 

70182587000 

64 

A$62 

70180037000 

65 

S$62 

70180038000 

66 

M$62 

70180039000 

67 

D$62 

70180040000 

68 

C$16 

70180059000 

69 

DBLE 

70180058000 

70 

CSQRT 

70182592000 

71 

CCOS 

70180066000 

72 

CSIN 

70182595000 

73 

CLOG 

70182591000 

74 

CEXP 

70182593000 

75 

CABS 

70182596000 

76 

E$5 1 

70182594000 

77 

A$52 

70180041000 

78 

S$52 

70180042000 

79 

M$52 

70180043000 

80 

D$52 

70180044000 

81 

A$55 

70182544000 

82 

S$55 

70180093000 

83 

M$55 

70182545000 

84 

D$55 

70180034000 

85 

CONJG 

70182598000 

86 

C$25 

70180068000 

87 

CMPLX 

70182597000 

88 

N$55 

END 

70180069000 
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MAGNETIC TAPE 70182806000 - LTCS2S 
(LIBRARY SOURCES CODED IN DAP) 


NUMBER 

NAME 

DOC. NO. 

1 

DMAX1 

70182585000 

2 

DMIN1 

70182586000 

3 

DINT 

70180850000 

4 

Z$80 

70180851000 

5 

A$81 

70180852000 

6 

C$61 

70182554000 

7 

A$66 

70180853000 

8 

A$66XRA 

70180979000 

9 

H$66 

70180855000 

10 

C$26 

70180857000 

11 

H$55 

70180860000 

12 

MAXO 

70182548000 

13 

MAXI 

70182549000 

14 

MINO 

70180649000 

15 

MINI 

70182551000 

16 

TANH 

70182565000 

17 

SQRT 

70182560000 

18 

SQRTX 

70180681000 

19 

SIN, COS 

70182563000 

20 

ATAN 

70182564000 

21 

E$2 1 

70182562000 

22 

E$22 

70180045000 

23 

ALOG 

70182559000 

24 

ALOGX 

70180682000 

25 

EXP 

70182561000 

26 

E$1 1 

70182547000 

27 

E$1 IX 

70180684000 

28 

ABS 

70182570000 

29 

C$62 

70180884000 

30 

AMOD 

70182572000 

31 

L$66 

70180854000 

32 

AINT 

70182571000 

33 

N$66 

70180856000 

34 

DIM 

70182573000 

35 

SIGN 

70182574000 

36 

AIM AG 

70180858000 

37 

L$55 

70180859000 

38 

IFIX 

70182553000 

39 

FLOAT 

70180062000 

40 

C$12 

70182575000 

41 

C$21 

70182558000 

42 

LOC 

70181962000 

43 

C$81 

70180882000 

44 

ISTORE 

70181982000 

45 

N$33 

70180090000 

46 

IFETCH 

70181983000 

47 

IABS 

70182552000 

48 

F$OE 

70181984000 

49 

MOD 

70182555000 

50 

F$TR-RA 

70180827000 

51 

SUB$ 

70185150000 
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FILE NUMBER 

NAME 

DOC. NO. 

5 2 

FbGA 

/U16DI0 iUUU 

5 3 

FSGC 

70] 85 15 2000 

54 

IDIM 

70182556000 

55 

A $2 2 

70182536000 

56 

M$22 

70182537000 

57 

A$22X 1 1 

70181805000 

58 

M$22X 1 1 

70181806000 

59 

D$22X 1 1 

70181804000 

60 

I SIGN 

70182557000 

61 

L$22 

70182534000 

62 

H$22 

70182535000 

63 

N$22 

70180097000 

64 

SLITE 

70182599000 

65 

M$ 1 1 

70180035000 

66 

D$1 1 

70182546000 

67 

M$ 1 IX 

70180685000 

68 

D$1 IX 

70180686000 

69 

OVERFL 

70180894000 

70 

F$AT 

70180071000 

71 

L$33 

70180065000 

72 

F$WN 

70180089000 

7? 

FtPN 

70180088000 

74 

F$R1 

70182610000 

75 

F$W 1 

70182611000 

76 

F$R2 

70182612000 

77 

FSW 2 

70182613000 

78 

F$R3 

70182614000 

79 

F$W 3 

70181667000 

80 

F$W 4 

70182616000 

81 

F$R5 - 9 

70180306000 

82 

F$F5-9 

70180310000 

83 

F$W5-9 

70180307000 

84 

F$10 

70182618000 

85 

016CHAIN 

70180659000 

86 

ARG$ 

70180072000 

87 

F $D 5 - 9 

70180308000 

88 

F$B5- 9 

70180309000 

89 

F$ER-RA 

70181068000 

90 

AC 1 

END 

70180717000 

MAGNETIC TAPE 70182803541 

- LTCM1S 

(LIBRARY OBJEC 

TS - SOFTWAR 

E VERSION) 


This magnetic tape consists of the concatenation of the individual objects of the listed 
programs. They have been translated by the FORTRAN Translator Mod 1 if FORTRAN coded 
and/or assembled by the DAP- 16 Mod 2 Assembler. 


FILE NUMBER 

NAME 

DOC. NO. 

1 

STMEAN 

70181386000 

2 

STGEOM 

70181387000 

3 

STCORR 

70181388000 

4 

STSPER 

70181389000 

5 

STCRMT 

70181390000 

6 

STMEDT 

70181391000 

7 

STCH12 

70181392000 
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FILE NUMBER 

NAME 

DOC. NO. 

8 

STBNPB 

70181394000 

9 

STLNRG 

70181395000 

10 

STCHIS 

70181393000 

11 

STPLRG 

70181396000 

12 

STANV1 

70181397000 

13 

STANV2 

70181398000 

14 

STANVR 

70181399000 

15 

STANVL 

70181400000 

16 

STANVG 

70181401000 

17 

STANYB 

70181402000 

18 

STANVY 

70181403000 

19 

DEFOAD 

70181405000 

20 

DEFOMA 

70181404000 

21 

DEFOHA 

70181406000 

22 

DEFORK 

70181407000 

23 

DESOAD 

70181408000 

24 

DESOMA 

70181409000 

25 

DESOHA 

70181410000 

26 

DESORK 

70181411000 

27 

PLYMUL 

70181416000 

28 

PLYDIV 

70181417000 

29 

PLYINT 

70181420000 

30 

PLYIRT 

70181421000 

31 

OLYEVL 

70181422000 

32 

PLYDIF 

70181423000 

33 

NACPLY 

70181429000 

34 

NAAITK 

70181418000 

35 

NALAGR 

70181419000 

36 

NABAIR 

70181424000 

37 

MATTRS 

70181412000 

38 

MATMUL 

70181431000 

39 

MARITH 

70181414000 

40 

MATEIG 

70181415000 

41 

MATINV 

70181427000 

42 

NAREGU 

70181425000 

43 

NAMULL 

70181426000 

44 

GSEID 

70181428000 

45 

SORT 

70181430000 

46 

SORT 2 

70181431000 

47 

CVPOLR 

70181432000 

48 

E$62 

70180053000 

49 

E$6l 

70180052000 

50 

E$26 

70182582000 

51 

E$66 

70180054000 

52 

DSQRT 

70182580000 

53 

DCOS 

70180055000 

54 

DSIN 

70182583000 

55 

DEXP 

70182581000 

56 

DLOG10 

70180051000 

57 

DLOG 

70182579000 

58 

DLOG2 

70182914000 

59 

DATAN2 

70180056000 

60 

DATAN 

70182584000 

61 

DMOD 

70182588000 

62 

DSIGN 

70182589000 

63 

DABS 

70182587000 

64 

A$62 

70180037000 
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NAME 

DOC. MO. 

65 

SS62 

70180036000 

66 

MS62 

70180039000 

67 

D$62 

70180040000 

68 

C$16 

70180059000 

69 

DELE 

70180058000 

70 

CSQRT 

70182592000 

71 

CCOS 

70180066000 

72 

CSIN 

70182595000 

73 

CLOG 

70182591000 

74 

CEXP 

70182593000 

75 

CABS 

70182596000 

76 

E$5 1 

70182594000 

77 

A$52 

70180041000 

78 

S$52 

70180042000 

79 

M$52 

70180043000 

80 

D$52 

70180044000 

81 

A$55 

70182544000 

82 

S$55 

70180093000 

83 

M$55 

70182545000 

84 

D$55 

70180034000 

85 

CONJG 

70182598000 

86 

C$25 

70180068000 

87 

CMPLX 

70182597000 

88 

N$55 

70180069000 

89 

DMAX1 

70182585000 

90 

DMIN1 

70182586000 

91 

DINT 

70180850000 

92 

Z$80 

70180851000 

93 

A$81 

70180852000 

94 

C$61 

70182554000 

95 

A$66 

70180853000 

96 

H$66 

70180855000 

97 

C$26 

70180857000 

98 

H$55 

70180860000 

99 

MAXO 

70182548000 

100 

MAXI 

70182549000 

101 

MINO 

70180649000 

102 

MINI 

70182551000 

103 

TANH 

70182565000 

104 

SQRT 

70182560000 

105 

SIN, COS 

70182563000 

106 

ATAN 

70182564000 

107 

E$2 1 

70182562000 

108 

E$22 

70180045000 

109 

ALOG 

70182559000 

110 

EXP 

70182561000 

111 

E$ll 

70182547000 

112 

ABS 

70182570000 

113 

C$62 

70180884000 

114 

AMOD 

70182572000 

115 

L$66 

70180854000 

116 

AINT 

70182571000 

117 

N$66 

70180856000 

118 

DIM 

70182573000 

119 

SIGN 

70182574000 

120 

AIMAG 

70180858000 
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FILE NUMBER 

NAME 

DOC. NO. 

121 

L$55 

70180859000 

122 

IF IX 

70182553000 

123 

FLOAT 

70180062000 

124 

C$12 

70182575000 

125 

C$21 

70182558000 

126 

LOC 

70181962000 

127 

C$81 

70180882000 

128 

ISTORE 

70181982000 

129 

N$33 

70180090000 

130 

IFETCH 

70181983000 

131 

IABS 

70182552000 

132 

F$DE 

70181984000 

133 

MOD 

70182555000 

134 

F$TR-RA 

70180827000 

135 

SUB$ 

70185150000 

136 

F$GA 

70185151000 

137 

F$GC 

70185152000 

138 

ID IM 

70182556000 

139 

A$22 

70182536000 

140 

M$22 

70182537000 

141 

ISIGN 

70182557000 

142 

L$22 

70182534000 

143 

H$22 

70182535000 

144 

N$22 

70180097000 

145 

SLITE 

70182599000 

146 

M$ 1 1 

70180035000 

147 

D$ 1 1 

70182546000 

148 

OVERFL 

70180894000 

149 

F$AT 

70180071000 

150 

L$33 

70180065000 

151 

F$WN 

70180089000 

152 

F$RN 

70180088000 

153 

F$R 1 

70182610000 

154 

F$W 1 

70182611000 

155 

F$R2 

70182612000 

156 

F$W2 

70182613000 

157 

F$R3 

70182614000 

158 

F$W 3 

7Q181667000 

159 

F$W4 

70182616000 

160 

F$R5-9 

70180306000 

161 

F$F5-9 

70180310000 

162 

F $ W 5-9 

70180307000 

163 

F$10 

70182618000 

164 

016CHAIN 

70180659000 

165 

ARG$ 

70180072000 

166 

F$D5-9 

70180308000 

167 

F$B5- 9 

70180309000 

168 

F$ER-RA 

70181068000 

169 

AC 1 

70180717000 

170 

SQRX1 

70188775000 

171 

COSX1 

70188781000 

172 

SINX1 

70188777000 

173 

ANT XI 

70188779000 

174 

LGEX1 

70188814000 

175 

LG2X1 

70188784000 

176 

EXEX1 

70188786000 

177 

EX2X1 

70188782000 
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FILE NUMBER 

NAME 

DOC. NO. 

178 

DSQRX1 

70 188788000 

179 

DCOSX1 

70188792000 

180 

DSINX1 

70188790000 

181 

DAT NX 1 

70188793000 

182 

DLGEX1 

70188801000 

183 

DLG2X1 

70188795000 

184 

DEXEX1 

70188799000 

185 

DEX2X1 

70188797000 

186 

DMPY 

70188808000 

187 

MPY 

70188811000 

188 

DIV 

70188810000 

189 

DADD 

70188812000 

190 

DSUB 

70188813000 

191 

ROND 

70188805000 

192 

RODD 

70188804000 

193 

TWOS 

70188803000 


END 
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1-47 
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FORTRAN Library 
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MAGNETIC TAPE 70 1 8280454 1 - LTCM 1H 
(LIBRARY OBJECTS - HARDWARE VERSION) 


This magnetic tape consists of the concatenation of the individual objects of the listed 
programs. They have been translated by the FORTRAN Translator Mod 1 if FORTRAN codf 
and/or assembled by the DAP-16 Mod 2 Assembler. 


NUMBER 

NAME 

DOC. NO. 

1 

STMEAN 

70181386000 

2 

STGEON 

70181387000 

3 

STCORR 

70181388000 

4 

STSPER 

70181389000 

5 

STCRMT 

70181390000 

6 

STMEDT 

70181391000 

7 

STCHI2 

70181392000 

8 

STBNPB 

70181394000 

9 

STLNRG 

70181395000 

10 

STCHIS 

70181393000 

11 

STPLRG 

70181396000 

12 

STANV 1 

70181397000 

13 

STANV2 

70181398000 

14 

STANVR 

70181399000 

15 

STANVL 

70181400000 

16 

STANVG 

70181401000 

17 

STANVB 

70181402000 

18 

STANVY 

70181403000 

19 

DEFOAD 

70181405000 

20 

DEFOMA 

70181404000 

21 

DEFOHA 

70181406000 

22 

DEFORK 

70181407000 

23 

DESOAD 

70181408000 
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FILE NUMBER 

NAME 

DOC. NO. 

24 

DESONA 

70181409000 

25 

DESOHA 

70181410000 

26 

DESORK 

70181411000 

27 

PYLMUL, 

70181416000 

28 

PLYDIV 

70181417000 

29 

PLYINT 

70181420000 

30 

PLYIRT 

70181421000 

31 

PLYEVL 

70181422000 

32 

PLYDIF 

70181423000 

33 

NACPLY 

70181429000 

34 

NAAITK 

70181418000 

35 

NALAGR 

70181419000 

36 

NAB AIR 

70181424000 

37 

NATTRS 

70181412000 

38 

NATMUL 

70181413000 

39 

MARITH 

70181414000 

40 

NATEIG 

70181415000 

41 

NATINY 

70181427000 

42 

NAREGU 

70181425000 

43 

NAMULL 

70181426000 

44 

GSEID 

70181428000 

45 

SORT 

70181430000 

46 

SORT 2 

70181431000 

47 

CVPOLR 

70181432000 

48 

E$62 

70180053000 

49 

E$6l 

70180052000 

50 

E$26 

70182582000 

51 

E$66 

70180054000 

52 

DSQRT 

70182580000 

53 

DCOS 

70180055000 

54 

DSIN 

70182583000 

55 

DEXP 

70182581000 

56 

DLOG10 

70180051000 

57 

DLOG 

70182579000 

58 

DLOG2 

70182914000 

59 

DATAN2 

70180056000 

60 

DATAN 

70182584000 

61 

DMOD 

70182588000 

62 

DSIGN 

70182589000 

63 

DABS 

70182587000 

64 

A$62 

70180037000 

65 

S$62 

70180038000 

66 

M$62 

70180039000 

67 

D$62 

70180040000 

68 

C$16 

70180059000 

69 

DBLE 

70180058000 

70 

CSQRT 

70182592000 

71 

CCOS 

70180066000 

72 

CSIN 

70182595000 

73 

CLOG 

70182591000 

74 

CEXP 

70182593000 

75 

CABS 

70182596000 

76 

E$5 1 

70182594000 

77 

A$52 

70180041000 

78 

S$52 

70180042000 

79 

M$52 

70180043000 

80 

D$52 

70180044000 
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NUMBER 

MAMK 

DOC. NO. 

81 

A $5 5 

70182544000 

82 

S$55 

70180003000 

83 

M$55 

70182545000 

84 

D$55 

70180034000 

85 

CONJG 

70182598000 

86 

C$25 

70180068000 

87 

CMPLX 

70182597000 

88 

N$55 

70180069000 

89 

DMAX1 

70182585000 

90 

DMIN1 

70182586000 

91 

DINT 

70180850000 

92 

Z$80 

70180851000 

93 

A$81 

70180852000 

94 

C$61 

70182554000 

95 

A$66X1\A 

70180979000 

96 

H$66 

70180855000 

97 

C$26 

70180857000 

98 

H$55 

70180860000 

99 

MAXO 

70182548000 

100 

MAXI 

70182549000 

101 

Ml NO 

70180640000 

102 

MINI 

70182551000 

103 

TANH 

70182565000 

104 

SQRTX 

70180681000 

105 

SIN, COS 

70182560000 

106 

ATAN 

70182564000 

107 

E$21 

70182562000 

108 

E$22 

70180045000 

109 

ALOGX 

70180682000 

110 

EXP 

70182561000 

111 

E$1 IX 

70180684000 

112 

ABS 

70182570000 

113 

C$62 

70180884000 

114 

AMOD 

70182572000 

115 

L$66 

70180854000 

116 

AINT 

70182571000 

117 

N$66 

70180856000 

118 

DIM 

70182573000 

119 

SIGN 

70182874000 

120 

AIMAG 

70180858000 

121 

L$55 

70180859000 

122 

IF IX 

70182553000 

123 

FLOAT 

70180062000 

124 

C$12 

70182575000 

125 

C$21 

70182558000 

126 

LOC 

70181962000 

127 

C$81 

70180882000 

128 

ISTORE 

70181982000 

129 

N$33 

70180090000 

130 

IFETCH 

70181983000 

131 

IABS 

70182552000 

132 

F$OE 

70181984000 

133 

MOD 

70182555000 

134 

F$TR-RA 

70180827000 

135 

SUB$ 

70185150000 

136 

F$GA 

70185151000 

137 

F$GC 

70185152000 
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FILE NUMBER 

NAME 

DOC. NO. 

138 

IDIM 

70182556000 

139 

A$22X1 1 

70181805000 

140 

M$22X 1 1 

70181806000 

141 

D$22X 1 1 

70181804000 

142 

ISIGN 

70182557000 

143 

L$22 

70182534000 

144 

H$22 

70182535000 

145 

N$22 

70180097000 

146 

SLITE 

70182599000 

147 

M$1 IX 

70180685000 

148 

D$1 IX 

70180686000 

149 

OVERFL 

70180894000 

150 

F$AT 

70180071000 

151 

L$33 

70180065000 

152 

F$WN 

70180089000 

153 

F$RN 

70180088000 

154 

F$R1 

70182610000 

155 

F$W 1 

70182611000 

156 

F$R2 

70182612000 

157 

F$W2 

70182613000 

158 

F$R3 

70182614000 

159 

F$W 3 

70181667000 

160 

F$W 4 

70182616000 

161 

F$R5-9 

70180306000 

162 

F$F5-9 

70180310000 

163 

F$W5-9 

70180307000 

164 

F$I0 

70182618000 

165 

016CHAIN 

70180659000 

166 

ARG$ 

70180072000 

167 

F $D 5-9 

70180308000 

168 

F$B5-9 

70180309000 

169 

F$ER-RA 

70181068000 

170 

AC 1 

70180717000 

171 

SQRX2 

70188776000 

172 

COSX2 

70180761000 

173 

SINX2 

70188778000 

174 

ATNX2 

70188780000 

175 

LGEX2 

70188815000 

176 

LG2X2 

70188785000 

177 

EXEX2 

70188787000 

178 

EX2X2 

70188783000 

179 

DSQRX2 

70188789000 

180 

DCOSX2 

70180762000 

181 

DSINX2 

70188791000 

182 

DATNX2 

70188794000 

183 

DLGEX2 

70188802000 

184 

DLG2X2 

70188796000 

185 

DEXEX2 

70188800000 

186 

DEX2X2 

70188798000 

187 

DM.PYH 

70188809000 

188 

DADD 

70188812000 

189 

DSUB 

70188813000 

190 

ROND 

70188805000 

191 

RODD 

70188804000 

192 

TWOS 

END 

70188803000 
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F iles 

1-47 

Statistical Library 

F iles 

48-170 

FORTRAN Library 

Files 

171-192 

Fixed Point Math Library 


PAPER TAPE 70181876000 - LTCF1 (Paper Tape 1 of 6) 

This paper tape consists of the concatenation of the individual objects of the listed 

the FORTRAN Translator MOD 1 and assembled by 


programs. They have been translated by 
the DAP- 16 Mod 2 Assembler. 

FILE NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


NAME 

DOC. NO. 

E$62 

70180053000 

E$6l 

70180052000 

E$26 

70182582000 

E$66 

70180054000 

DSQRT 

70182580000 

DCOS 

70180055000 

DSIN 

70182583000 

DEXP 

70182581000 

DLOG10 

70180051000 

DLOG 

70182579000 

DLOG2 

70182914000 

DATAN2 

70180056000 

DATAN 

70182584000 

DMOD 

70182588000 

DSIGN 

70182589000 

DABS 

70182587000 

A$62 

70180037000 

S$62 

70180038000 

M$62 

70180039000 

D$62 

70180040000 

C$16 

70180059000 

DBLE 

70180058000 

END 



PAPER TAPE 70181877000 - LTCF2 (Paper Tape 2 of 6) 

This paper tape consists of the concatenation of the individual objects of the listed 
program. They have been translated by the FORTRAN Translator MOD 1 and assembled by 
the DAP-16 Mod 2 Assembler. 


FILE NUMBER 

NAME 

DOC. NO. 

1 

CSQRT 

70182592000 

2 

CCOS 

70180066000 

3 

CSIN 

70182595000 

4 

CLOG 

70182591000 

5 

CEXP 

70182593000 

6 

CABS 

70182596000 

7 

E$5 1 

70182594000 

8 

A$52 

70180041000 

9 

S$52 

70180042000 

10 

M$52 

70180043000 
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FILE NUMBER 

NAME 


DOC. NO. 

11 

D$52 


70180044000 

12 

A$55 


70182544000 

13 

S$55 


70180093000 

14 

M$55 


70182545000 

15 

D$55 


70180034000 

1 6 

CONJG 


70182598000 

17 

C$25 


70180068000 

18 

CMPLX 


70182597000 

19 

N$55 


70180069000 


END 



PAPER TAPE 

70181882000 - LTCF3S 

(Paper 

Tape 3 of 6) 

the Software Version of Paper Tape 3 



FILE NUMBER 

NAME 


DOC. NO. 

1 

DMAX1 


70182585000 

2 

DMIN1 


70182586000 

3 

DINT 


70180850000 

4 

Z$80 


70180851000 

5 

A$81 


70180852000 

6 

C$61 


70182554000 

7 

A$66 


70180853000 

8 

H$66 


70180855000 

9 

C$26 


70180857000 

10 

H$55 


70180860000 

11 

MAX0 


70182548000 

12 

MAXI 


70182549000 

13 

MI NO 


70180649000 

14 

MINI 


70182551000 

15 

TANH 


70182565000 

16 

SQRT 


70182560000 

17 

SIN, COS 


70182563000 

18 

ATAN 


70182564000 

19 

E$2 1 


70182562000 

20 

E$22 


70180045000 

21 

ALOG 


70182559000 

22 

EXP 


70182561000 

23 

E$1 1 


70182547000 


END 



PAPER TAPE 

70181878000 - LTCF3H 

(Paper 

Tape 3 of 6) 


This is the Hardware Version of Paper Tape 3 


FILE NUMBER 

NAME 

DOC. NO. 

1 

DMAX1 

70182585000 

2 

DMIN1 

70182586000 

3 

DINT 

70180850000 

4 

Z$80 

70180851000 

5 

A$81 

70180852000 

6 

C$61 

70182554000 

7 

A$66XRA 

70180979000 

8 

H$66 

70180855000 

9 

C$26 

70180857000 
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FILE NUMBER 

NAME 

DOC. NO. 

10 

H$55 

70180860000 

1 1 

MAXO 

70182548000 

12 

MAXI 

70182549000 

13 

MINO 

70180649000 

14 

MINI 

70182551000 

15 

TANH 

70182565000 

16 

SQRTX 

70180681000 

17 

SIN, COS 

70182563000 

18 

ATAN 

70182564000 

19 

E$2 1 

70182562000 

20 

E$22 

70180045000 

21 

ALOGX 

70180682000 

22 

EXP 

70182561000 

23 

E$1 IX 

70180684000 


END 


PAPER TAPE 

70181879000 - LTCF4 (Paper 

Tape 4 of 6) 

FILE NUMBER 

NAME 

DOC. NO. 

1 

ABS 

70182570000 

2 

/-■•'Wo 

0^)06 

70 180884000 

3 

AMOD 

70182572000 

4 

L$66 

70180854000 

5 

AINT 

70182571000 

6 

N$66 

70180856000 

7 

DIM 

70182573000 

8 

SIGN 

70182574000 

9 

AIMAG 

70180858000 

10 

L$55 

70180859000 

11 

IFIX 

70182553000 

12 

FLOAT 

70180062000 

13 

C$12 

70182575000 

14 

C$21 

70182558000 

15 

LOC 

70181962000 

16 

C$81 

70180882000 

17 

ISTORE 

70181982000 

18 

N$33 

70180090000 

19 

IF ETCH 

70181983000 

20 

LABS 

70182552000 

21 

F$OE 

70181984000 

22 

MOD 

70182555000 

23 

F$TR-RA 

70180827000 

24 

SUB$ 

70185150000 

25 

F$GA 

70185151000 

26 

F$GC 

70185152000 


END 


PAPER TAPE 

70181883000 - LTCF5S (Paper Tape 5 of 6) 


This is the Software Version of Tape 5 


FILE NUMBER 

NAME 

DOC. NO. 

1 

IDEM 

70182556000 

2 

A$22 

70182536000 

3 

M$22 

70182537000 

4 

ISIGN 

70182557000 
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FILE NUMBER 

NAME 

DOC. NO. 

5 

L$22 

70182534000 

6 

H$22 

70182535000 

7 

N$22 

70180097000 

8 

SLITE 

70182599000 

9 

M$1 1 

70180035000 

10 

D$1 1 

70182546000 

11 

OVERFL 

70180894000 

12 

F$AT 

70180071000 

13 

L$33 

70180065000 

14 

F$WN 

70180089000 

15 

F$RN 

70180088000 

16 

F$R1 

70182610000 

17 

F$W 1 

70182611000 

18 

F$R2 

70182612000 

19 

F$W2 

70182613000 

20 

F$R3 

70182614000 

21 

F$W 3 

70181667000 

22 

F$W4 

70182616000 

23 

F$R5-9 

70180306000 

24 

F$F5- 9 

70180310000 

25 

F $W 5-9 

70180307000 

26 

F$IO 

70182618000 


END 



PAPER TAPE 70181880000 - LTCF5H (Paper Tape 5 of 6) 
This is the Hardware Version of Tape 5 


NUMBER 

NAME 

DOC. NO. 

1 

IDIM 

70182556000 

2 

A$22X 1 1 

70181805000 

3 

M$22X1 1 

70181806000 

4 

D$22X1 1 

70181804000 

5 

ISIGN 

70182557000 

6 

L$22 

70182534000 

7 

H$22 

70182535000 

8 

N$22 

70180097000 

9 

SLITE 

70182599000 

10 

M$1 IX 

70180685000 

11 

D$1 IX 

70180686000 

12 

OVERFL 

70180894000 

13 

F$AT 

70180071000 

14 

L$33 

70180065000 

15 

F$WN 

70180089000 

16 

F$RN 

70180088000 

17 

F$R1 

70182610000 

18 

F$W 1 

70182611000 

19 

F$R2 

70182612000 

20 

F$W2 

70182613000 

21 

F$R3 

70182614000 

22 

F$W 3 

70181667000 

23 

F$W4 

70182616000 

24 

F$R5-9 

70180306000 

25 

F$F5- 9 

70180310000 

26 

F $W 5-9 

70180307000 

27 

F$IO 

70182618000 


END 
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PAPER TAPE 70181881000 - LTCF6 (Paper Tape 6 of 6) 


FILE NUMBER NAME 


1 

2 

3 

4 

5 

6 


016CHAIN 

ARG$ 

F$D5-9 

F$B5-9 

F$ER-RA 

AC 1 

END 


DOC. NO. 

70180659000 

70180072000 

70180308000 

70180309000 

70181068000 

70180717000 
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APPENDIX B 


MATHEMATICAL ROUTINES 


F unction Routine 

Complex 


Absolute value 

CABS 

Add 

A$55 

Add real argument 

A$52 

Conjugate 

CONJG 

Convert imaginary part to real 

AIMAG 

Cosine 

CCOS 

Divide 

D$55 

Divide by real argument 

D$52 

Exponential, base e 

CEXP 

Load 

L$55 

Load real part 

REAL 

Logarithm, base e 

CLOG 

Multiply 

M$55 

Multiply by real argument 

M$52 

Negate 

N$55 

Raise to integer power 

E$5 1 

Sine 

CSIN 

Square root 

CSQRT 

Store (hold) 

H$55 

Subtract 

S$55 

Subtract single-precision argument 

S$52 

Double -Precis ion 


Absolute value 

DABS 

Add 

A$66 

Add single-precision argument 

A$62 

Add integer to exponent 

A$81 

Arctangent, principal value 

DATAN 
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F unction 


Routine 


Double - Precision 


Arctangent, X/Y' 

DATAN2 

Clear (zero) exponent 

Z$80 

Convert exponent to integer 

C$81 

Convert to integer 

C$61 

Convert to single -precis ion 

C$62 

Cosine 

DCOS 

Divide 

D$66 

Divide by real argument 

D$62 

Exponential, base e 

DEXP 

Load 

L$66 

Logarithm, base e 

DLOG 

Logarithm, base 2 

DLOG2 

Logarithm, base 10 

DLOGIO 


DM AX 1 

Minimum value 

DMIN1 

Multiply 

M$66 

Multiply by real argument 

M$62 

Negate 

N$66 

Raise to double-precision power 

E$66 

Raise to integer power 

E$6l 

Raise to real power 

E$62 

Remainder 

DMOD 

Sine 

DSIN 

Square root 

DSQRT 

Store (hold) 

H$66 

Subtract 

S$66 

Subtract real argument 

S$62 

Transfer sign of second argument to first 

DSIGN 

Truncate fractional bits 

DINT 

Truncate fractional bits and convert 
to integer 

IDINT 


Real 


Absolute value 

ABS 

Add 

A$22 

Arctangent, principal value 

ATAN 

Arctangent, X/Y 

ATAN2 

Convert pair to complex 

CMP LX 
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F unction 


Routine 


Real 


Convert to complex format 

C$25 

Convert to double-precision 

C$26 

Convert to integer 

C$21 

Divide 

D$22 

Exponential, base e 

EXP 

Hyperbolic tangent 

TANH 

Load 

L$22 

Logarithm, base e 

ALOG 

Logarithm, base 10 

ALOGIO 

Maximum integer value 

MAXI 

Maximum value 

AMAX1 

Minimum integer value 

MINI 

Minimum value 

AMIN1 

Multiply 

M$22 

Positive difference 

DIM 

Raise to double-precision power 

E$26 

Raise to integer power 

E$2 1 

Raise to real power 

E$22 

Remainder 

AMOD 

Sine, cosine 

SIN, COS 

Square root 

SORT 

Store (hold) 

H$22 

Subtract 

S$22 

Transfer sign of second argument to first 

SIGN 

T runcate fractional bits 

AINT 

Truncate fractional bits and convert 
to integer 

IF IX, INT 

TWOs complement 

N$22 


Integer 


Absolute value 

IABS 

Convert to double-precision 

C$16 

Convert to real (FORTRAN-generated) 

FLOAT 

Convert to real 

C$12 

Divide 

D$1 1 

Maximum value 

AMAXO 


B-3 


AM 74 



Function 


Routine 


Maximum integer value MAXO 

Minimum value AMINO 

Minimum integer value MINO 

Multiply M$ll 

Positive difference IDIM 

Raise to integer power E$ll 

Remainder MOD 

Transfer sign of second argument to first ISIGN 


Logical 


Complement 
OR with A-register 


N$33 

L$33 
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APPENDIX C 


SUBROUTINE FUNCTIONS 


INTRINSIC AND EXTERNAL FUNCTIONS 


Mathematical and Trigonometric Functions 



Argument 

Result 



Name 

Data Type 

Data Type 

Function 


SIN 

R 

R 

Sine (radians) 


DSIN 

D 

D 



CSIN 

C 

C 



COS 

R 

R 

Cosine (radians) 


DCOS 

D 

D 



CCOS 

C 

C 



AT AN 

R 

R 

Arctangent (radians) 


DATAN 

D 

D 



ATAN2 

R 

R 



DATAN2 

D 

D 



TANH 

R 

R 

Hyperbolic tangent (radians) 


SQRT 

R 

R 

Square root 


DSQRT 

D 

D 



CSQRT 

C 

C 



EXP 

R 

R 

Exponential 


DEXP 

D 

D 



CEXP 

C 

C 



ALOG 

R 

R 

Natural logarithm 


DLOG 

D 

D 


CLOG 

C 

C 



ALOG 10 

R 

R 

Common logarithm 


DLOG2 

D 

D 


DLOGIO 

D 

D 



ABS 

R 

R 

Absolute value 


IABS 

I 

I 



DABS 

D 

D 



CABS 

C 

R 



AMOD 

R 

R 

Remainder 


MOD 

I 

I 



DMOD 

D 

D 



AINT 

R 

R 

Truncate fractional bits 


DINT 

D 

I 



IDINT 

D 

I 



IFIX 

R 

I 



INT 

R 

I 
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Argument Result 


Name 

Data Type 

Data Type 

Function 

AMAXO 

I 

R 

Choose largest argument 

AMAX1 

R 

R 


MAXO 

I 

I 


MAXI 

R 

I 


DMAX1 

D 

D 


AMINO 

1 

R 

Choose smallest argument 

AMIN 1 

R 

R 


DM1N1 

D 

D 


MIN 10 

I 

I 


MINI 

R 

I 


FLOAT 

I 

R 

Change data type or argument 

AIMAG 

C 

R 


DBLE 

R 

D 


CMPLX 

C 

R 


REAL 

C 

R 


SNGL 

R 

D 


SIGN 

R 

R 

Value of first argument, sign of 

DSIGN 

D 

D 

second 

1SIGN 

I 

I 


DIM 

R 

R 

Positive difference 

IDIM 

I 

I 


CON JG 

C 

C 

Complex conjugate 


necial Subroutines for FORTRAN Use 


Argument Result 

Name Data Type Data Type 

IFETCH(I) 

ISTORE(I, J) 

LOC 

OVERFL 

SLITE 

SLITET 

SSWTCH 


Function 

Get contents of location I 
Store value of J in location I 
Find address of argument 
Check for error condition 
Set and reset sense lights or 
switches 


COMPILER SUPPORT SUBROUTINES 


Conversion Routines 



Argument 

Result 


Name 

Data Type 

Data Type 

Function 

C$12 

I 

R 

Convert integer to real 

C$16 

I 

D 

Convert integer to double -precision 

C$21 

R 

I 

Convert real to integer 

C$25 

R 

C 

Convert real to complex 

C$26 

R 

D 

Convert real to double -precision 

C$61 

D 

I 

Convert double -precision to integer 

C$62 

D 

R 

Convert double -precision to real 

C$81 

D 

D 

Convert exponent of double - 
precision number to integer 


C-2 


AM74 



Arithmetic Routines 


Name 

Function 

Name 

Function 

A$22 

R = R+R 

E$62 

D = D**R 

A$52 

C = C+R 

E$66 

D = D**D 

A$55 

C = C+C 

M$1 1 

I = 1*1 

A$62 

D = D+R 

M$22 

R = R*R 

A$66 

D = D+D 

M$52 

C = C*R 

A$81 

D = D*(2**I) 

M$55 

C = C*C 

D$1 1 

I = I/I 

M$62 

D = D*R 

D$22 

R = R/R 

M$66 

D = D*D 

D$52 

C = C/R 

N$22 

R = -R 

D$55 

C = C/C 

N$33 

L = -L 

D$62 

D = D/R 

N$55 

C = -C 

D$66 

D = D/D 

N$66 

D = -D 

E$1 1 

I = 1**1 

S$22 

R = R-R 

E$2 1 

R = R**I 

S$52 

C = C-R 

E$22 

R = R**R 

S$55 

C = C-C 

E$26 

D = R**D 

S$62 

D = D-R 

E$5 1 

C = C**I 

S$66 

D = D-D 

E$6 1 

D = D**I 

Z$80 

Replace binary exponent 
with zero 


Miscellaneous Routines 

Name Function 

AC 1 Pseudo accumulators 

ARG$ Convert indirect address to direct address 

H$22 Store real number in memory 

H$55 Store complex number in memory 

H$66 Store double -precision number in memory 

L$22 Load real number into A- and B -registers 

L$55 Load complex number into complex accumulator 

L$66 Load double -precision number into double -precision accumulator 

L$33 INCLUSIVE OR with A-register 

SUB$ Calculate address of array element 


C-3 


AM74 



APPENDIX D 


LIBRARY INDEX 


Primary 

Name 

Entry 

Points 

Subroutines 

Called 

A$22 

A$22 

ARG$ 


S$22 

N$22 

F$ER 

A$22X 

A$22 

N$22 


S$22 

F$ER 

A$52 

A$52 

F$AT 

H$55 

L$22 

A$22 

H$22 

L$55 

A$55 

A$55 

F$AT 

H$55 

SUB$ 

L$22 

A$22 

H$22 

L$55 

A$62 

A$62 

F$AT 

H$66 

DBLE 

A$66 

A$66 

A$66 

N$66 


S$66 

F$ER 


M$66 

H$66 


D$66 

L$66 

ARG$ 

AC 1 

AC2 

AC 3 

A$66X 

A$66 

N$66 


A$66X 

F$ER 


S$66 

H$66 


S$66X 

L$66 


M$66 

ARG$ 


M$66X 

AC 1 


D$66 

AC 2 


D$66X 

AC 3 


Approx. 

Number of Storage Tape 

References (Words ) Number 


1 150 5S 

1 
1 

1 140 5H 

1 

1 20 2 


1 

1 

1 

1 

1 

1 60 2 

1 
4 
2 
2 
2 
1 

1 20 1 

1 
1 
1 

11 580 3S 

3 
1 
1 
1 
1 
1 
1 

11 530 3H 

3 
1 
1 
1 
1 
1 
1 


D- 1 
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Primary 

Name 

Entry 

Points 

Subroutines 

Called 

Number of 
References 

Approx. 
Storage 
(Words ) 

T ape 
Number 

A$81 

A$81 

N$22 

2 

70 

3 



F$ER 

1 





AC 1 

1 





AC2 

1 



ABS 

ABS 

L$22 

1 

10 

4 



N$22 

1 



AC 1 

AC 1 



5 

6 


AC2 






AC 3 






AC4 






AC 5 





AIMAG 

AIMAG 

L$55 

1 

10 

4 



L$22 

1 





AC 3 

1 



AINT 

AINT 

L$22 

1 

30 

4 



N$22 

2 





A$22 

1 





S$22 

1 



A LOG 

ALOGIO 

ARG$ 

1 

120 

3S 


A LOG 

C$12 

1 





H$22 

5 





L$22 

3 





A$22 

6 





S$22 

2 





D$22 

1 





M$22 

7 





F$ER 

1 



ALOGX 

ALOGIO 

ARG$ 

1 

180 

3H 


ALOG 

C$12 

1 




ALOGX 

A$22 

4 





M$22 

4 





S$22 

1 





F$ER 

1 



ALOGIO 

See ALOG 

or ALOGX 




AMAXO 

See MAXO 





AMAX1 

See MAXI 





AMINO 

See MINO 





AMIN1 

See MINI 





A MOD 

AMOD 

L$22 

1 

30 

4 



D$22 

1 





AINT 

1 





M$22 

1 





N$22 

1 





A$22 

1 



ARG$ 

ARG$ 



20 

6 


D-2 


AM74 



Primary- 

Name 

Entry 

Points 

Subroutines 

Called 

Number of 
References 

Approx. 

Sto rage 
(Words 

T ape 
Number 

ATAN 

ATAN2 

ARG$ 

3 

340 

3 


ATAN 

D$22 

6 





N$22 

7 





M$22 

5 





A$22 

11 





S$22 

2 



ATAN2 

See ATAN 




C$12 

C$12 

A$22 

1 

30 

4 



N$22 

1 



C$16 

C$16 

C$12 

1 

5 

1 



C$26 

1 



C$21 

C$21 

N$22 

1 

30 

4 



A$22 

1 





F$ER 

1 



C$25 

C$25 

H$22 

1 

20 

2 



CMPLX 

1 



C$26 

C$26 

AC 1 

1 

10 

3 



AC2 

1 





AC 3 

1 



C$61 

C$61 

C$62 

1 

4 

3 



C$21 

1 



C$62 

C$62 

L$22 

1 

20 

4 


SNGL 

N$66 

1 





N$22 

1 





L$66 

1 





AC 1 

1 





AC2 

1 



C$81 

C$81 

AC 1 

1 

10 

4 

CABS 

CABS 

F$AT 

1 

40 

2 



SUB$ 

2 





L$22 

2 





M$22 

2 





H$22 

2 





A$22 

1 





SQRT 

1 



CCOS 

CCOS 

F$AT 

1 

40 

2 



L$55 

1 





A$55 

1 





H$55 

1 





CSIN 

1 



CEXP 

CEXP 

F$AT 

1 

60 

2 



SUB$ 

7 





EXP 

1 





H$22 

3 





COS 

1 





M$22 

2 





SIN 

1 





L$55 

1 




D-3 
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Primary 

Name 

Entry 

Points 

Subroutines 

Called 

Number of 
References 

Approx. 
Storage 
(Words iq) 

tape 

Number 

CLOG 

CLOG 

F$AT 

1 

90 

2 



SUB$ 

6 





L$22 

3 





M$22 

3 





H$22 

5 





A$22 

1 





ALOG 

1 





ATAN2 

1 





L$55 

1 



CMP LX 

CMPLX 

F$AT 

I 

40 

2 



SUB$ 

2 





L$22 

2 





H$22 

2 





L$55 

1 



CONJG 

CONJG 

F$AT 

1 

40 

2 



SUB$ 

4 





L$22 

2 





H$22 

2 





N$22 

1 



COS 

See SIN 

L$55 




CSIN 

CSIN 

F$AT 

1 

90 

2 



SUB$ 

5 





EXP 

1 





H$22 

6 





L$22 

3 





D$22 

1 





A$22 

1 





SIN 

1 





M$22 

4 





S$22 

1 





COS 

1 





L$55 

1 



CSQRT 

CSQRT 

F$AT 

1 

90 

2 



SUB$ 

7 





CABS 

1 





H$22 

8 





ABS 

1 





A$22 

1 





M$22 

2 





SQRT 

1 





L$22 

6 





D$22 

1 





L$55 

1 



D$1 1 

D$1 1 

ARG$ 

1 

80 

5S 



F$ER 

1 



D$1 IX 

D$1 1 

ARG$ 

1 

40 

5H 


D$1 IX 

F$ER 

1 



D$22 

See M$22 





D$22X 

D$22 

N$22 

3 

110 

5H 



F$ER 

2 






D-4 


AM74 



Primary- 

Name 

D$52 


D$55 


D$62 


D$66 

D$66X 

DABS 

DATAN 


DATAN2 


DBLE 


Entry Subroutines Number of 

Points Called References 


D$52 F$AT 

H$55 
SUB$ 
L$22 
D$22 
H$22 
L$55 

D$55 F$AT 

H$55 
SUB$ 
L$22 
M$22 
H$22 
A$22 
D$22 
S$22 
N$22 
L$55 

D$62 F$AT 

H$66 
DBLE 
L$66 
D$66 

See A$66 

See A$66X 

DABS F$AT 

L$66 
N$66 


1 

1 

2 

2 

2 

2 

1 

1 

1 

12 

8 

6 

8 

2 

2 

1 

1 

1 

1 

2 

1 

1 

1 


1 

1 

1 


DATAN F$AT 1 

DABS 1 

H$66 9 

C$81 1 

L$66 1 3 

A$66 10 

N$66 3 

D$66 2 

M$66 9 

DATAN2 F$AT 1 

L$66 9 

H$66 3 

F$ER 1 

D$66 1 

DATAN 1 

S$66 1 

A$66 1 

DBLE F$AT 1 

L$22 1 

C$26 1 


Approx. 

Storage Tape 

(Words ) Number 

50 2 


140 2 


20 1 


10 1 


180 1 


70 1 


20 1 


D-5 
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Primary 

Name 

Entry 

Points 

Subroutines 

Called 

Number of 
References 

Approx. 
Storage 
(Words 1Q ) 

T ape 
Number 

DCOS 

DCOS 

F$AT 

1 

20 

1 



LS66 

1 





A$66 

1 





H$66 

1 





DSIN 

1 



DEXP 

DEXP 

F$AT 

1 

160 

1 



L$66 

12 





M$66 

8 





H$66 

11 





C$61 

1 





C$16 

1 





N$66 

1 





A$66 

8 





S$66 

3 





D$66 

1 





A$8 1 

1 



DIM 

DIM 

L$22 

1 

20 

4 



S$22 

1 



DINT 

DINT 

L$66 

1 

20 

3 



N$66 

2 





A$66 

1 





S$66 

1 





AC 1 

1 



DIV $ 

See M$22 





DLOG 

DLOG 

F$AT 

1 

10 

1 



DLOG2 

1 





M$66 

1 



DLOG 2 

DLOG2 

F$AT 

1 

100 

1 



L$66 

5 





F$ER 

1 





C$81 

1 





C$16 

1 





H$66 

6 





Z$80 

1 





A$66 

6 





S$66 

2 





D$66 

1 





M$66 

6 



DLOG 10 

DLOG 10 

F$AT 

1 

10 

1 



DLOG2 

1 





M$66 

1 



DMAX1 

DMAX1 

L$66 

3 

40 

3 



H$66 

2 





S$66 

1 



DMIN1 

DMIN1 

L$66 

3 

40 

3 



H$66 

2 





S$66 

1 




D-6 
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Primary- 

Name 

Entry 

Points 

Subroutines 

Called 

Number of 
References 

Approx. 
Storage 
(Words 1Q ) 

Tape 

Number 

DMOD 

DMOD 

F$AT 

1 

20 

1 



L$66 

1 





D$66 

1 





H$66 

1 





DINT 

1 





M$66 

1 





S$66 

1 





N$66 

1 



DSIGN 

DSIGN 

F$AT 

1 

20 

1 



L$66 

3 





N$66 

1 



DSIN 

DSIN 

F$AT 

1 

130 

1 



DABS 

1 





M$66 

9 






H$66 

C$61 

C$16 

N$66 

A$66 

MOD 

L$66 

S$66 

DSQRT 

DSQRT 

F$AT 

L$66 

C$62 

H$22 

SQRT 

C$26 

H$66 

D$66 

A$66 

A$81 

E$1 1 

E$1 1 

ARG$ 
M$1 1 
F$ER 

E$1 IX 

E$1 1 

arg$ 


E$1 IX 

F$ER 

E$2 1 

E$21 

ARG$ 

M$22 

D$22 

E$22 

E$22 

ARG$ 

ALOG 

M$22 

EXP 

E$26 

E$26 

F$AT 

C$26 

H$66 

DLOG 

M$66 

DEXP 


5 

1 

1 

3 

7 
1 

8 
2 

1 40 1 

2 
1 
1 
1 
1 
1 
1 
1 
1 

1 100 3S 

2 

1 

1 110 3H 

1 

1 50 3 

1 
1 

1 30 3 

1 
1 
1 

1 30 1 

1 

2 
1 
1 
1 


D-7 
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Primary 

Name 

Entry 

Points 

Subroutines 

Called 

Number of 
.References 

Approx. 
Storage 
fvV urds q) 

T ape 

in uiiiioe i 

E$5 1 

E$5 1 

F$AT 

1 

60 

2 



H$55 

3 





LABS 

1 





L$55 

4 





M$55 

1 





D$55 

1 



E$6l 

E$6l 

F$AT 

1 

70 

1 



H$66 

5 





L$66 

5 





D$66 

1 





D$1 1 

2 





M$1 1 

1 





M$66 

2 



E$62 

E$62 

F$AT 

1 

30 

1 



H$66 

2 





DLOG 

1 





M$62 

1 





DEXP 

1 



E$66 

E$66 

F$AT 

1 

30 

1 



H$66 

2 





DLOG 

1 





M$66 

1 





DEXP 

1 



EXP 

EXP 

ARG$ 

1 

230 

3 



N$22 

2 





M$22 

6 





S$22 

3 





A$22 

2 





D$22 

2 





F$ER 

1 



F$AR 

See F$IO 




5 

F$AT 

F$AT 



58 

5 

F$B5-9 

F$B5 

C$MR 

6 

26 

6 


F$B6 





- 

F$B7 






F$B8 






F$B9 






F$BN 





F$CB 

See F$IO 




5 

F $D 5-9 

F$D2 

o$ps 

2 

34 

6 


F$D5 

0$ME 

6 




F$D6 

F$D7 

F$D8 

F$D9 

F$DN 


D-8 
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Primary- 

Name 

Entry 

Points 

Subroutines 

Called 

Number o£ 
References 

Approx. 
Storage 
(Words jq) 

Tape 

Number 

F$ER 

F$ER 

ACS 

2 

37 

6 


F$HT 





F$F5-9 

F$F5 

C$BR 

1 

41 

5 


F$F 6 

F$ER 

1 




F$F7 






F$F8 






F$F9 






F$FN 





F$GA 

F$GA 

F$ER 

1 

18 

4 

F$GC 

F$GC 



14 

4 

F$HT 

See F$ER 




6 

F$IO 

F$IO 

F$ER 

2 

1356 

5 


F$CB 

F$L1 

F$L2 

F$L3 

F$L5 

F$L6 

F$AR 


F$R1 

F$R1 

F$IO 

1 

21 

5 



I$AA 

1 





I$AB 

1 



F$R2 

F$R2 

F$IO 

1 

21 

5 



I$PA 

1 





I$PB 

1 



F$R3 

F$R3 

F$IO 

1 

21 

5 



I$CA 

1 





I$CB 

1 



F$R5-9 

F$R5 

F$IO 

1 

80 

5 


F$R6 

I$MA 

1 




F$R7 

I$MC 

1 




F$R8 






F$R9 





F$Rn 

F$Rn 

F$R1 

1 

45 

5 



F$R2 

1 





F$R3 

1 





F$R5 

1 





F$R6 

1 





F$R7 

1 





F$R8 

1 





F$R9 

1 



F$TR 

F$TR 

F$W 1 

5 

198 

4 



F$AR 

4 





F$CB 

1 





F$L6 

1 





AC 1 

2 





AC 2 

1 





AC 3 

1 




D-9 



r* ima ry 

Name 

Flntry 

Points 

Subroutines 

Called 

Number of 
References 

Approx. 
Storage 
(Words ) 

T ape 
Numbe r 

F $W 1 

F$W 1 

F$IO 

1 

94 

5 



0$AP 

1 





0$AC 

1 





0$AF 

1 





0$AB 

1 



F$W2 

F$W2 

F$IO 

1 

80 

5 



0$PF 

1 





o$pp 

1 





o$pc 

1 





0$PB 

1 



F$W 3 

F$W 3 

F$IO 

1 

39 

5 



0$CH 

1 





0$CB 

1 



F$W 4 

F$W 4 

F$IO 

1 

36 

5 



0$LF 

3 





0$LP 

1 





0$L0 

1 



F $W 5-9 

F$W 5 

F>1G 

1 

57 

5 


F$W 6 

0$MC 

1 




F $W 7 

C$8T06 

1 




F$W8 

0$MA 

1 




F$W 9 





F$Wn 

F$Wn 

F$W 1 

1 

41 

5 



F$W2 

1 





F$W 3 

1 





F$W 4 

1 





F$W 5 

1 





F$W 6 

1 





F$W7 

1 





F$W8 

1 





F$W 9 

1 



FLOAT 

FLOAT 

C$12 

1 

10 

4 

H$22 

H$22 

ARG$ 

1 

10 

5 

H$55 

H$55 

ARG$ 

1 

20 

3 



AC 1 

1 





AC2 

1 





AC 3 

1 





AC4 

1 



H$66 

H$66 

ARG$ 

1 

20 

3 



AC 1 

1 





AC2 

1 





AC 3 

1 



IABS 

IABS 



10 

4 

IDIM 

IDIM 



20 

5 

IDINT 

See IF IX 





IFETCH 

IFETCH 

ARG$ 

1 

10 

4 


D - 10 


AM 74 



Approx. 

Primary Entry Subroutines Number of Storage Tape 

Name Points Called References (Words „) Number 


IF IX 

ID INT 

L$22 


INT 

IF IX 

C$21 

INT 

See IF IX 


ISIGN 

ISIGN 


ISTORE 

ISTORE 

F$AT 

L$22 

REAL 

L$22 

ARG$ 

L$33 

L$33 


L$55 

L$55 

ARG$ 
AC 1 

AC2 

AC 3 

AC4 

L$66 

L$66 

ARG$ 
AC 1 

AC2 

AC3 

LOC 

LOC 


M$1 1 

M$1 1 

ARG$ 

F$ER 

M$1 IX 

M$1 1 

ARG$ 


M$1 IX 

F$ER 

M$22 

M$22 

N$22 


D$22 

ARG$ 


DIV $ 

F$ER 

M$22X 

M$22 

F$ER 

M$52 

M$52 

F$AT 


H$55 

SUB$ 

L$22 

M$22 

H$22 

L$55 

M$55 M$55 F$AT 

H$55 

SUB$ 

L$22 

M$22 

H$22 

S$22 

N$22 

A$22 

L$55 


1 10 4 

1 

20 5 

1 10 4 

1 10 5 

10 5 

1 20 4 


1 

1 

1 

1 


1 20 4 

1 
1 
1 

10 4 

1 110 5S 

1 

1 50 5H 

1 

5 330 5S 

2 
3 

1 130 5H 

1 50 2 


1 

2 

2 

2 

2 

1 

1 110 2 
1 

10 

4 

4 

4 

1 

1 

1 

1 


D-ll 
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Approx. 


P rimary 

Entry 

Subroutines 

Number of 

Storage 

T ape 

Name 

Points 

Called 

References 

(Words 1Q ) 

Number 

M$62 

M$62 

F$AT 

1 

20 

1 



H$66 

1 





DBLE 

1 





M$66 

1 



M$66 

See A$66 




35 

M$66X 

See A$66X 




34 

MAXO 

A MAXO 

MAXO 

FLOAT 

1 

40 

3 

MAXI 

AMAX1 

L$22 

2 

50 

3 


MAXI 

H$22 

2 





S$22 

1 





IFIX 

1 



MINO 

AMINO 

MINO 

FLOAT 

1 

30 

3 

MINI 

AMIN1 

L$22 

2 

50 

3 


MINI 

H$22 

2 





S$22 

1 





IFIX 

1 



MOD 

MOD 

D$1 1 

1 

20 

4 



M$1 1 

1 



N$22 

N$22 



10 

5 

N$33 

N$33 



10 

4 

N$55 

N$55 

H$55 

1 

30 

2 



SUB$ 

2 





L$22 

2 





N$22 

2 





H$22 

2 





L$55 

1 



N$66 

N$66 

AC 1 

1 

30 

4 



AC2 

1 





AC 3 

1 



OVERFL 

OVERFL 

AC 5 

1 

20 

5 

REAL 

See L$22 




5 

S$22 

See A$22 




5S 

S$22X 

See A$22X 





S$52 

S$52 

F$AT 

1 

30 

2 



H$55 

1 





L$22 

1 





S$22 

1 





H$22 

1 





L$55 

1 



S$55 

S$55 

F$AT 

1 

40 

2 



H$55 

1 





SUB$ 

4 





L$22 

2 





S$22 

2 
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3 rimary 

Name 

Entry- 

Points 

Subroutines 

Called 

Number of 

References 

Approx. 
Storage 
(Words i o ) 

T ape 
Number 



N$22 

2 





H$22 

2 





L$55 

1 



S$62 

S$62 

F$AT 

1 

20 

1 



H$66 

1 





DBLE 

1 





S$66 

1 





N$66 

1 



S$66 

See A$66 




3S 

S$66X 

See A$66X 




3H 

SIGN 

SIGN 

L$22 

2 

20 

4 



N$22 

1 



SIN 

COS 

ARG$ 

1 

190 

3 


SIN 

N$22 

2 





M$22 

7 





S$22 

1 


- 



A$22 

4 



SLITE 

SLITE 

ARG$ 

3 

70 

5 


SLITET 

L$33 

1 




SSWTCH 





SLITET 

See SLITE 





SNGL 

See C$62 





SQRT 

SQRT 

ARG$ 

1 

70 

3S 



DIV $ 

1 





D$22 

1 





A$22 






F$ER 




SQRTX 

SQRT 

ARG$ 

1 

80 

3H 


SQRTX 

D$22 

1 





A$22 

1 





F$ER 

1 



SSWTCH 

See SLITE 





SUB$ 

SUB$ 

M$1 1 

3 

130 

4 



F$ER 

1 



TANH 

TANH 

L$22 

1 

60 

3 



EXP 

1 





A$22 

2 





H$22 

1 





D$22 

1 



Z$80 

Z$80 

AC 1 

1 

20 

3 


D - 1 3 


A 



APPENDIX E 


ERROR MESSAGES 


Error 


Message 

Condition 

Subroutine 

AD 

Over /underflow in double -precision 

A$66, S$66, A$66X, S$66X 

BF 

End-of-file mark encountered while unit 
backspacing a record. 

F$F5 -9 

DL 

Negative or zero argument 

DLOG, DLOGIO, DLOG2 

DT 

Both arguments are zero 

DATAN2 

DZ 

Division by zero 

D$22, D$22X 

EQ 

Exponential overflow adding integer to 
double -precision exponent 

A$81 

EX 

Exponential overflow during exponentiation 

EXP 

FE 

Format error 

F$IO 

GO 

Incorrect control variable in a GO TO 

statement 

F$GA 

II 

First argument zero, second argument 
negative 

I > 2 and J > 15, or 

I < -2 and J > 1 5 

E$ 1 1 , E$1 IX 

IM 

Over /underflow during integer 
mul ti pli c atio n 

M$1 1 , M$1 IX 

IN 

Input error 

F$AR 

IZ 

Integer division by zero or -32, 768/ -1 

D$1 1 , D$1 IX 

LG 

Log of negative or zero argument 

ALOG, ALOGIO, ALOGX 

MD 

Double -precision multiplication or 
division over /underflow 

D$66, M$66, D$66X, M$66X 

PZ 

Double-precision division by zero 

D$66, D$66X 

RI 

Integer too large when converted from 
real to integer 

C$21 

SA 

Arithmetic overflow (result > 2**127) 

A$22, A$22X 

SD 

Divisor unnormalized 

D$22 

SM 

Arithmetic overflow during multiplication 
or division 

M$22, M$22X, D$22X 

SQ 

Negative argument 

SQRT, SQRTX 
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